Draft version July 2, 2012 

Preprint typeset using I4TgX style emulateapj v. 08/22/09 



SPITZER IRAC IDENTIFICATION OF HERSCHEL-ATLAS SPIRE SOURCES 

Sam Kim', Julie L. Wardlow', Asantha Cooray', S. Fleuren^, W. Sutherland^, A. A. Khostovan', R. Auld', M. Baes"*, 
R. S. Bussmann\ S. Buttiglione*, a. Cava^, D. Clements**, A. Dariush**, G. DeZotti*', L. Dunne'", S. Dye", S. Eales\ 
J. Fritz", R. Hopwood*, E. Ibar'', R. Ivison'^, M. Jarvis''' '^, S. Maddox'", M. J. Michalowski'\ E. Pascale\ M. Pohlen\ 
E. Rigby"', D. ScoTT'^ D. J. B. Smith'*, R Temi", and P. van der Were'* 

Draft version July 2, 2012 

ABSTRACT 

We use Spitzer-IRAC data to identify near-infrared counterparts to submilHrneter galaxies detected with 
He rschel-SPlRE at 250 /im in the Herschel Astrophysical Terahertz Large Area Survey (H- ATLAS). Using a 
likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159 in the 
survey area. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 
3.6 and 4.5 /im color-magnitude space, and we use this property to identify further 23 counterparts to 13 SPIRE 
sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated 
with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification 
rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 
have both 3.6 and 4.5 /im coverage. Of these, 43 have [3.6] — [4.5] > indicating that they are likely to be 
at z > 1.4. Thus, ^ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We 
discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, 
dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have lO^'' to 
10^^ Lq total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a 
factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples. 
Subject headings: galaxies: high-redshift - infrared: galaxies - galaxies:starburst 



1. INTRODUCTION 

The extragalactic background at far-infrared (IR) and sub- 
millimeter (sub-mm) wavelengths is well-constrained from 
total intensity measurements (Puget et al. 1996; Fixsen et 
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al. 1998; Dwek et al. 1998). However, the properties of the 
discrete galaxies that make up this background are still largely 
unknown. These sub-mm galaxies are expected to capture the 
dusty star-formation out to redshifts of 4 and beyond and are 
now understood to be an integral component of galaxy forma- 
tion and evolution (Hughes et al. 1998; Eales et al. 1999; Blain 
et al. 2002; Chapman et al. 2003; Austermann et al. 2010). 

A challenge for studies of sub-mm galaxies is that wide- 
field surveys at these wavelengths have poor spatial reso- 
lution, caused by the limited apertures of single dish sub- 
mm survey telescopes (Small, Ivison & Blain 1997; Scott 
et al. 2002; Coppin et al. 2006). Typical resolutions are of 
the order of tens of arcseconds, making the identification 
of multiwavelength counterparts challenging, particularly in 
the optical where sub-mm galaxies are usually faint due to 
high dust extinction. Furthermore, the surface density of 
faint optical galaxies is such that several potential counter- 
parts may be situated within the positional error of each sub- 
mm source. However, the identification of counterparts to 
sub-mm galaxies is critical for both photometric and spectro- 
scopic studies, which yield redshifts, spectral energy distribu- 
tions (SEDs) and morphological information. The measure- 
ments of these quantities for statistically significant samples 
of sub-mm galaxies are required to compare their properties 
with theoretical predictions, and to fully understand the role 
of the sub-mm bright phase in galaxy evolution. 

The most reliable way to pinpoint the positions of sub-mm 
galaxies is through high-resolution sub-mm interferometry, 
which directly traces the dust emission at wavelengths com- 
parable to the original selection function (e.g. Downes et al. 
1999; Gear et al. 2000; lono et al. 2006; Younger et al. 2007). 
The sensitivity of early generations of sub-mm interferomet- 
ric instrumentation limited anlayses to samples of no more 
than a few sources (e.g. Downes et al. 1999; Gear et al. 2000; 
lono et al. 2006; Younger et al. 2007, 2008, 2009; Ivison 
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et al. 2008; Cowie et al. 2009; Aravena et al. 2010a; Hat- 
sukade et al. 2010; Tamura et al. 2010; Chen et al. 2011 
), and although rapid progress is being made (e.g. Smolcic 
2012a,b) such observations are still unfeasible for large sur- 
veys of sources. In addition, extended sources (e.g. Matsuda 
et al. 2007), or sources that are a blend of multiple compo- 
nents (e.g. Wang et al. 201 1; Smolcic 2012b), can be difficult 
to detect. 

An alternative technique to identify sub-mm galaxies is to 
take advantage of the far-lR-radio correlation (e.g. Condon et 
al. 1998; Garrett 2002) and the relatively low surface density 
of radio sources by searching for counterparts in deep inter- 
ferometric radio data (e.g. Ivison et al. 1998, 2002; Borys 
et al. 2004; Dye et al. 2009; Dunlop et al. 2010; Biggs et 
al. 2011). Similarly, the dusty, active sub-mm galaxies are 
typically bright at mid-lR wavelengths, a property that has 
also been utilized to identify counterparts (e.g. Ivison et al. 
2004; Pope et al. 2006; Dye et al. 2008, 2009; Clements et al. 
2008; Dunlop et al. 2010; Biggs et al. 201 1). Deep radio and 
mid-lR (typically Spitzer MIPS 24 yumbut also Spitzer IRAC 
3.6 /^tm; Biggs et al. 2011) searches typically identify coun- 
terparts to 60-80% of sub-mm galaxies, and nearest-neighbor 
positional matching is then used to investigate the properties 
of these sources at other wavelengths. The remaining 20- 
40% of sources are thought to be at high-redshifts (z > 3) or 
be dominated by cold dust (To < 25 K at z ~ 2; Chapman et 
al. 2005). 

The Herschel Astrophysical Terahertz Large Area Survey 
(H- ATLAS; Bales et al. 2010) is the largest open-time key 
project being undertaken by the Herschel Space Observa- 
tory (Pilbratt et al. 2010)'''. The total planned survey area 
is 550 deg^ within which it is expected to detect > 300,000 
bright sub-mm galaxies. During the Science Demonstration 
Phase (SDP), H-ATLAS observed 14.4 deg^ in the Galaxy 
And Mass Assembly (GAMA; Driver et al. 201 1) 9 hour field, 
to 5(7 depths of 35 to 132 mJy in five bands from 100 to 
500 ^m, using PACS (100 and l6Q^im; Poglitsch et al. 2010) 
and SPIRE (250, 350 and 5Q0^lm\ Griffin et al. 2010) in par- 
allel mode. The lack of deep radio and 24 /im data across the 
field restricts counterpart identification to optical (e.g. Sloan 
Digital Sky Survey [SDSS]; York et al. 2000) or near-lR data 
(e.g. VISTA Kilo-Degree Infrared Galaxy Survey [VIKING]; 
Sutherland et al. in prep, UKIRT Infrared Deep Sky Survey 
[UKIDSS]; Lawrence et al. 2007). 

Smith et al. (2011) identified SDSS and GAMA galaxies 
as the counterparts to 37% of the SPIRE sources in the H- 
ATLAS SDP field. In the wider area of the whole GAMA- 
9 hour field Fleuren et al. (2012) increased the fraction of 
SPIRE sources with reliable identifications to 51% by using 
A's-band imaging from the VIKING survey. In the GAMA- 
15 hour field 50% of the SPIRE sources have reliable coun- 
terparts identified in Wide-Field Infrared Survey Explorer 
(WISE) data at 3.4 ^m (Bond et al. 2012). 

In this paper we cross-identify H-ATLAS SPIRE sources 
with Spitzer Infrared Array Camera (IRAC; Fazio et al. 
2004a) galaxies, which are selected from 3.6 and 4.5-/im 
observations of 0.4 deg^ of the H-ATLAS SDP field. The 
IRAC data are advantageous for counterpart identification be- 
cause of their mid-lR wavelengths and depth (~3 magnitudes 
deeper than WISE). Thus, counterpart identification is less bi- 

" Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with important 
participation from NASA. 



ased towards galaxies with relatively low redshifts or dust ob- 
scurations than the existing analyses. 

The paper is organized as follows: the data analysis and 
sample selection are presented in Section |2] In Section [3] 
we describe the counterpart identification method and in Sec- 
tion 2] we discuss our results, including the identification rate 
and the colors and properties of SPIRE counterparts. Our con- 
clusions are presented in Section |5] A table of the identified 
counterparts, including magnitudes and fluxes is presented in 
the Appendix. We use J2000 coordinates and ACDM cosmol- 
ogy with n„^= 0.27, riA= 0.73 and = 70kms-iMpc^i 
throughout. All photometry is on the AB magnitude system 
where 1 /iJy= 23.9 mag; IRAC 3.6 and 4.5 /im AB magni- 
tudes are designated [3.6] and [4.5], respectively. 

2. SAMPLE SELECTION AND DATA ANALYSIS 

The H-ATLAS fields were chosen to minimize contamina- 
tion from Galactic cirrus and to maximize the overlap with ex- 
isting and planned wide-area imaging and spectroscopic sur- 
veys. The SDP field, which we study here, spans ^ 14 deg^ 
in the GAMA-9 hour field and has existing SDSS, VIKING, 
and UKIDSS data. It is observed with PACS at 100 and 
160 Aim and SPIRE at 250, 350 and 500 ^im. The H-ATLAS 
PACS and SPIRE map-making processes are described in Ibar 
et al. (2010) and Pascale et al. (2011), respectively, and de- 
tails of the source extraction procedures are given in Rigby 
et al. (2011). A public catalog of SDP sources is available 
from the H-ATLAS website-". The data reach 5<j point-source 
depths of 132, 126, 32, 36 and 45 mJy, with beam sizes of 8.7, 
13.1, 18.1, 25.2 and 36.9" (FWHM) at 100, 160, 250, 350, 
and 500 ^m, respectively. Although PACS 100 and 160/im 
data have significantly higher angular resolution than SPIRE 
data we do not use PACS data for counterpart identification 
because only five of the 159 SPIRE sources studied here are 
detected by PACS. In all five cases the position of the iden- 
tified IRAC counterpart (table [2]) is consistent with the PACS 
emission. Furthermore, the counterparts to all of these five 
sources are low-redshift late-type galaxies, which is consis- 
tent with expectations for galaxies that are detectable with 
PACS. 

Spitzer IRAC 3.6 and 4.5 staring mode observations of 
0.4 deg^, in three regions of the H-ATLAS SDP field were 
carried out during the warm mission (S/j/fzer program 548; PI: 
A. Cooray). Two of the areas targeted contain bright SPIRE 
SDP sources that are now known to be gravitationally lensed 
(ID81 and ID130; Negrello et al. 2010; Hopwood et al. 201 1) 
and the third area was chosen as a test field. In staring mode 
simultaneous observations at 3.6 and 4.5 /im are offset from 
each other by 6.8'. Due to the offset, 0.22 deg'^ of the tar- 
geted area has imaging data at both 3.6 and 4.5 /im, and the 
remaining 0.18 deg-^ is split between3.6 and 4.5 /im coverage. 
There are 159 SPIRE sources in the total IRAC footprint. Of 
these, 101 are observed at both 3.6 and 4.5 /im, and 30 (28) 
have only 3.6 /im (4.5 /im) data. The total exposure time is 
1080 seconds per pixel in each filter. Data reduction and mo- 
saicking is performed on the Corrected Basic Calibrated Data 
(cBCD) with MOPEX (MOsaicker and Point source Extrac- 
tor; Makovoz & Marleau 2005). The final mosaicked images 
have 0.6" X 0.6" pixels and the FWHM of the IRAC point 
spread function (PSF) at 3.6 and 4.5 /im is 1.66 and 1.72", 
respectively. 

http://www.h-atlas.org 
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Our source detection and extraction procedure is designed 
to simultaneously select sources at 3.6 and 4.5 /xm, and pro- 
vide equivalent photometry in both bands. The source detec- 
tion is performed on a mean combined 3.6 and 4.5 fim image. 
The 3.6 and 4.5 /.tm image combination is carried out with 
Montage^', which models background flux by co-adding 
and re-projecting background corrected images into a mosaic. 
Source detection is performed with S EXTRACTOR (Bertin & 
Arnouts 1996) and detected sources are required to comprise 
at least three contiguous pixels with fluxes at least 1.5tT above 
the local background. There are 27170 unique 3.6 and 4.5 /im 
sources that meet this requirement and that are detected at 
> 5ct. Photometry is measured in 3.8" diameter apertures 
with the APPHOT task in IRAF. The advantage of APPHOT is 
that the photometry is measured in fixed apertures at specified 
source positions. The measured aperture photometry is cor- 
rected to "total" fluxes, assuming point-source profiles, and 
using the calibration derived by the SWIRE team for IRAC 
data with multiplicative correction factors of 1.36 and 1.40 at 
3.6 and 4.5 /im, respectively (Surace et al. 2005). 

The completeness of the IRAC catalog is established by in- 
serting artificial galaxies with a maximum half-light radius of 
1.5" in the IRAC images. Our source extraction and photo- 
metric procedure is repeated and a source is considered to be 
recovered if the extracted position and magnitude are within 
1.5" and 5a, respectively, of the input values. The 50% com- 
pleteness limits are 22.5 mag (3.63 fiJy) and 22.2 mag (4.79 
/iJy) at 3.6 and 4.5 fim respectively. 

The far-IR emission detected by SPIRE traces emission 
from cold dust. Therefore, point sources in blank-field sur- 
veys, such as H-ATLAS, are primarily expected to be exter- 
nal galaxies and not stars in the Milky Way. Indeed, of the ~ 
6600 SPIRE sources in - 14 deg^ in the H-ATLAS SDP field 
only 78 are Galactic stars, and two are candidate debris disks 
(Thompson et al. 2010). Since stars dominate the IRAC cat- 
alog at the brightest fluxes and have magnitude distributions 
that are different to galaxies, it is necessary to remove them 
prior to counterpart identification (e.g. Smith et al. 201 1). Un- 
fortunately the angular resolution of IRAC (FWHM ~ 1 .7") is 
insufficient to morphologically distinguish between stars and 
high-redshift galaxies. Furthermore, while stars can be iden- 
tified in [3.6] — [4.5] color-magnitude space (e.g. Eisenhardt 
et al. 2004) we are limited to just single IRAC band data for 
^ 45% of the total IRAC catalog. 

Therefore, instead of using IRAC data alone to reliably 
identify stars, we make use of the stellar classification in 
SDSS. For this we match the IRAC catalog to SDSS DR7 
(Abazajian, et al. 2009) using positional information. ^ 40% 
of IRAC sources are matched with SDSS catalog within 
search radius 1 .5". Out of total 27 170 IRAC sources, there are 
4239 (16%) sources that are classified as stars by SDSS. Once 
these are removed, the final IRAC catalog has 22931 galaxies. 
This catalog is used for the SPIRE source identification and 
subsequent analysis. We note that Fazio et al (2004b) showed 
that faint stars are a only minor contribution to the IRAC pop- 
ulation, with < 10'^ mag^^ deg^^ for magnitudes fainter than 
16 mag at both 3.6 and 4.5 /im. Therefore, stars that are too 
faint to be detected as such by SDSS are not expected to affect 
our statistical analysis. 

In addition, SDSS QSO population statistics (Schneider et 
al. 2007) indicate that only 2-3 QSOs are expected in the 
IRAC area (< 1 in the counterpart search area around SPIRE 

^ ' http ://montage . ipac . caltech.edu 
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Fig. 1. — Galaxy number counts from H-ATLAS/Spfeer area compared 
with previous IRAC deep and wide fields: SDWFS (Ashby et al. 2009), 
First Look Survey (FLS; Fazio et al. 2004b) and EGS (Barmby et al. 2008). 
H-ATLAS galaxy counts are more consistent with SDWFS. The 50% IRAC 
completeness limits of our imaging at 22.5 and 22.2 mag at 3.6 and 4.5 fj,m, 
respectively, are indicated by vertical dash-dot-dot lines. 5ct detection limit 
of SDWFS with 3" aperture is 22.8 and 22.3 mag at 3.6 and 4.5 fim, respec- 
tively. 

centroids), and therefore the effect of QSOs on our statistics 
is also expected to be negligible. We conclude that the coun- 
terpart identification statistics are unlikely to be affected by 
unresolved source contamination in our input catalog. 

In figure [T] we compare the IRAC galaxy counts in H- 
ATLAS/S/j/fzer area with some of the previous IRAC deep 
and wide fields. H-ATLAS galaxy counts are more consis- 
tent with SDWFS (Ashby et al. 2009), although the variations 
between fields are likely to be due to cosmic aariance. 

3. SPITZER IRAC IDENTIFICATION OF SPIRE 
SOURCES 

The size of the SPIRE beam is 18.1" (FWHM) at 250 /im, 
and 25.2 and 36.9" at 350 and 500 /im, respectively. The 
spatial density of IRAC-bright galaxies is high enough that 
multiple sources may lie within the SPIRE beam. Typi- 
cally, this will include the true SPIRE counterpart and unas- 
sociated foreground and background galaxies. However, the 
surface density of 250- /im bright sources is sufficient that 
a single SPIRE source may be composed of emission from 
multiple galaxies. Therefore, choosing the nearest object as 
the counterpart of a given SPIRE source can be misleading. 
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Inst ead, we perform a likelihood ratio (LR) analysis (Sec- 
tion EU Sutherland & Saunders 1992; CiUegi et al. 2003; 
Brusa et al.2007; Smith et al.2011; Fleuren et al. 2012), 
which uses positional information and the magnitude distri- 
bution of counterparts and background sources, to identify 
SPIRE sources in the IRAC catalog. Furthermore, the LR 
analysis shows that galaxies associated with SPIRE sources 
occupy a distinct region of IR AC color-flux space (figure |4|l, 
which we use in Section 13.41 to identify additional counter- 
parts. 

We note that other statistical matching techniques have been 
used in astronomical implementations. For example, the cor- 
rected Poissonian probability (p-statistic; Downes et al. 1986) 
uses the surface density to calculate the probability of a source 
with the observed magnitude being detected within the radius 
investigated (e.g. Ivison et al. 2002, 2007; Pope et al. 2006; 
Chapin et al. 2009). However, this technique is most ap- 
propriate for catalogues in which the surface density is low 
(e.g. radio data) and favors counterparts that are brighter than 
the background population. Our IRAC data do not have low- 
surface density, and we do not wish to assume that the Her- 
schel counterparts are typically the brightest IRAC sources. 
Bayesian techniques (e.g. Budavari & Szalay 2008), which 
use a priori knowledge of the counterpart population to guide 
the identification process, could also be employed (e.g. Brand 
et al. 2006; Gorjian et al. 2008). However, we do not wish to 
bias our results by assuming a priori knowledge of the SPIRE 
population at IRAC wavelengths. The LR analysis is advanta- 
geous because it computes the intrinsic IRAC magnitude dis- 
tribution of the SPIRE sources from the data provided (see 
e.g. Fleuren et al. 2012). 

3.1. The Likelihood Ratio method 

The LR, L, is the ratio of the probability that the IRAC 
source is the correct SPIRE counterpart with the equivalent 
probability for an unassociated background source. Following 
Sutherland & Saunders (1992), L is calculated as: 



L = 



q{m)f{r) 
n{m) 



(1) 



where q{m) and n{m) are the normalized magnitude distri- 
butions of counterparts and background sources, respectively. 
The radial probability distribution of the separation between 
the SPIRE 250 fim and IRAC positions is denoted by /(r). 

We estimate n{m), the normalized magnitude distribution 
of background sources, by averaging the source counts in ar- 
eas of circle within 7.2" radius of 100 different realizations 
of 159 random positions. The radius is three times the posi- 
tional uncertainty of the SPIRE catalog. We use cTpos = 2.4 
± 0.9" (Smith et al. 2011) calculated in the H-ATLAS SDP 
field from the positional offsets between SDSS DR7 galaxies 
and 6621 SPIRE sources. The center of each random circle is 
required to be at least 15" away from the nearest SPIRE cen- 
troid to minimize possible contamination of real association. 
The quantity q{m), the normalized magnitude distribution of 
true IRAC counterparts to SPIRE sources cannot be directly 
derived. We empirically estimate q{ra) by first subtracting the 
magnitude distribution of background galaxies from the mag- 
nitude distribution of all IRAC sources within search radius. 

This results in an estimate of the magnitude distribution of 
sources that are in excess of the background, iVcxcoss(™), and 
assumed to be the true counterparts to SPIRE sources. Then 
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Fig. 2. — The magnitude distribution of IRAC galaxies at 3.6 (top) and 
4.5 /jm (bottom) witliin the 7.2" counterpart search radius around SPIRE 
sources. We calculate the distribution of background galaxies and subtract 
these to reveal an excess around SPIRE sources which peaks at 19.0 mag at 
both 3.6 and 4.5 fim . These are the magnitude distributions of detected IRAC 
galaxies that are associated with SPIRE sources. The magnitude distribution 
of background sources is obtained by averaging the counts within 7.2" of 
100 different realizations of 159 random positions that are each > 15" from 
the nearest SPIRE centroid. To reduce the Poisson fluctuations associated 
with the small sample of 159 SPIRE centroids, the magnitude distribution of 
IRAC counts within SPIRE error circles is averaged over two bins. The 50% 
Spitzer-IVJKC completeness limits of 22.5 and 22.2 mag at 3.6 and 4.5 /-tm, 
respectively, are indicated by vertical dash-dot-dot lines. All of the excess 
galaxies that are associated with SPIRE sources are brighter than these de- 
tection limits indicating that the IRAC data are sufficiently deep to detect 
counterparts to all typical SPIRE sources. 



qija) is calculated via: 

q{m) = 



,(m) 



Em ^cxcoss("l) 



X Qq, 



(2) 



where Qq is the fraction of true counterparts above the IRAC 
detection limit. The [3.6] and [4.5] distributions of n{m) and 



^excess (jTi) are shown in figure[2] 



The radial probability distribution, /(r), is given by 



exp 



(3) 



where r is the separation between the SPIRE 250 /^m and 
IRAC centroids, and dpos is the positional uncertainty of the 
SPIRE catalog. The median 250 /Ltm signal-to-noise ratio of 
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the SPIRE sources in our sample is 6.33 and only 43 (27%) 
have SNR> 8. Therefore, we ignore the improvement in o-pos 
for sources with high signal-to-noise ratio (Smith et al. 201 1). 

Each H-ATLAS source may have several potential coun- 
terparts. Therefore, we calculate the reliability, Rj, of each 
IRAC galaxy, j, in SPIRE counterpart search radius. Rj is 
the probability that the galaxy, j, is the correct IRAC counter- 
part to the SPIRE source. Following Sutherland & Saunders 
(1992): 



Rj 



E,i, + (i-Qo) 



(4) 



where i represents each IRAC source in the search radius. 

To accept a potential counterpart as reliable we require that 
Rj > 0.8, i.e. there is < 20% chance of a false association. 
This criterion was used by Smith et al. (201 1) for the counter- 
part search with SDSS data and our results on the identifica- 
tion rate with IRAC data can be directly compared to their re- 
sults based on optical imaging. We note that the choice of the 
exact limiting value of R does not strongly affect our conclu- 
sions. Of th e 123 SPIRE sources with R > 0.8 counterparts 
(section [33]) 107 (87%) have R > 0.9, and 101 (82%) have 
R > 0.95. 

For a catalog in which identified sources are required to 
have i? > 0.8 the overall false detection rate is 



7V(falsc) = XI " -^j) ■ 



(5) 



Ri>0.i 



3.2. Qo calculation 

In order to apply the LR method to the SPIRE and IRAC 
catalogs we must first compute Qo- Smith et al. (2011) cal- 
culated Qo from the total number of SDSS galaxies in the er- 
ror circles of SPIRE sources. However, this method requires 
some a priori chosen radius and is only valid if each SPIRE 
source is associated with only one SDSS galaxy. Instead, we 
make use of an alternative method, derived by Fleuren et al. 
(2012), which is radius independent and unaffected by the 
clustering of IRAC sources. Clustering can be a significant 
complication because SPIRE sources are expected to reside in 
overdense regions similar to environments of sub-mm galax- 
ies (e.g. Aravena et al. 2010b; Hickox et al. 2012). 

For the purposes of this discussion we define a SPIRE er- 
ror circle without a true IRAC counterpart as a "blank". Di- 
rectly observed blanks arise due to counterparts below the 
IRAC detection limit or beyond the search radius, r. How- 
ever, the total number of blanks must also be statistically cor- 
rected for contamination by unassociated foreground or back- 
ground galaxies. Thus, the true number of blanks (St) is 
equal to the number of observed SPIRE blanks (S{r)) and 
the number of SPIRE sources that are falsely matched with 
random background galaxies. We define Ro{r) as the num- 
ber of random positions that contain a IRAC galaxy while 
R{r) = N — Ro{r) is the number of random blanks. Then, 



St = S{r) + [St 



Ro(r) 

N ^ 



(6) 



with Ro{r) and R{r) calculated from 100 random realizations 
of N = 159 error circles, which are all located at least 15" 
away from SPIRE centroids. 

Qo is defined as the fraction of true counterparts above 
the IRAC detection limit. Thus, for an infinite search ra- 
dius, the rate of true blanks, St /N is simply equal to 1 — Qo. 
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Fig. 3. — Quantities employed in the calculation of Qo- the 
ratio of SPIRE blanks to random blanks, is fitted with 1 — QoP{^) in order 
to calculate Qo. The best fit values are Qo = 0.93 ± 0.09 and 0.92 ± 0.08 
at 3.6 and 4.5 /^m, respectively. We also show R(r)/N and S{r)/N. the 
rate of random positions with at least one IRAC source and the rate of SPIRE 
blanks, respectively. 



Fleuren et al. (2012) demonstrate that equation |6] can be rear- 
ranged to show that St/N ^ S{r) /R{r), and thus 1 - Qo = 
S{r)/R{r). 

However, there is also the possibiUty that the true coun- 
terpart is outside of the examined area and this must be ac- 
counted for in the Qo estimate. The probability that the real 
SPIRE source is outside the search radius can be derived ana- 
lytically from the normalized SPIRE source distribution, f{r) 
(equation[3]l. This probabiUty, F{r), is 



F{r) 



P{r') dr' = 1 



(7) 



where P{r) = 27rr/(r) (Fleuren et al. 2012). Assuming that 
the probability of a SPIRE source being a true blank (1 — Qo) 
and the probability that the detected counterpart is outside of 
the search radius (1 — i^(r)) are independent, it can be shown 
that the total probability that there is no counterpart is (see 
Fleuren et al. 2012) 



| = l-QoF(r). 



(8) 



However, we have already shown that St/N — S{r)/R{r). 
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Fig. 4. — Top: IRAC color-magnitude diagram. All LR counterparts with 
JJ > 0.8 and potential counterparts are highlighted. Potential counterparts ai'e 
IRAC galaxies within 7.2" of SPIRE sources that are not otherwise identified 
with the LR method. The LR counterparts typically have [3.6] — [4.5] > — 
0.4 and are brighter the background IRAC population. IRAC galaxies above 
the solid line and within 7.2" of SPIRE positions have a probability of be- 
ing random associations of < 20%. This discriminator is used to identify 
an addit ional 23 IRAC counterparts to 13 SPIRE sources ("color IDs"; sec- 
tion|33). Bottom: SPIRE color-flux diagram for all sources in the H-ATLAS 
SDP field. We highlight sources with identified IRAC counterparts using the 
LR and color-magnitude methods. Sources that have formal 500 /^m fluxes 
below the Icr detection limit are shown as lower limits and are excluded from 
the statistical analysis. A two dimensional KS test between the background 
population and the identified SPIRE sources yields the probabiHty that the 
two datasets are drawn from the same parent population of p = 0.172. Thus, 
we conclude that the identified galaxies ai'e not a strongly biased subsample 
of SPIRE sources, although there may be a small selection effect. Average 
error bars are shown at the top-left hand corner of both plots. 

Thus, using the observables S{r) and R{r) one can calcu- 
late Qq via 1 - QoF{r) S{r)/R{r). Figure [3] shows 
S{r)/R{r), R{r)/N, and S{r)/N, observed for search radii 
of one to 10". We calculate the constant, Qq, by fitting 
1 — QoF{R) to S{r)/R{r). minimization yields Qq ~ 
0.93 ± 0.09, 0.92 ± 0.08 at 3.6 and 4.5 ^m, respectively. 

3.3. LR counterparts 



We apply the LR technique outlined in section l3TTl and 13.21 
at 3.6 and 4.5 /im, to identify IRAC counterparts to 123 of 
159 H-ATLAS SPIRE sources (Tables [T] and |2]). q{m) and 
n{m) are calculated in a radius qf3(Tpos = 7.2" (although 
Qo is radius independent; section [3?2] | and only IRAC galax- 
ies within 7.2" of SPIRE sources are considered to be po- 



tential counterparts. This limit includes 99.7% of true coun- 
terparts, whilst minimizing the potential contamination from 
unassociated sources. The analysis is performed separately 
at 3.6 and 4.5 fim and any IRAC galaxy with R > 0.8 at ei- 
ther wavelength is considered to be a SPIRE counterpart. In 
addition, where both 3.6 and 4.5 /im data is available we com- 
bine the probabilities (R) from each wavelength, and include 
five counterparts that have a combined probability of being 
considered to be reliable identifications > 0.8, but have R < 
0.8 at 3.6 and 4.5 um individually. All 123 counterparts are 
presented in table|2] including 81 galaxies with both 3.6 and 
4.5 yum photometry. The false identification rate is 1.9% or 
approximately two sources (equation |5]i; in the case of coun- 
terparts that are identified at both 3.6 and 4.5 /im we use the 
highest of the two Rj in this calculation. 

The identification rate of 77% is significantly higher than 
that from optical (37% in SDSS; Smith et al. 201 1) or near-IR 
analyses (51% in VIKING; Fleuren et al. 2012) of H-ATLAS 
sources. This difference is likely to be primarily driven by 
the wavelengths of the study, which are less sensitive to K- 
correction and dust absorption than optical data. 

The counterpart identification rate is also significantly 
higher than the 50% obtained with WISE 3.4-/im data (Bond 
et al. 2012); this is true even if we only consider the 3.6 //m 
data, where 97 of the 129 SPIRE sources with 3.6 fim cover- 
age (75%) are identified. WISE is shallower than our IRAC 
observations - 19.7 mag at 3.4 /im compared to 22.5 mag 
at 3.6 /im in IRAC. However, 79% of the 3.6-/im identified 
counterparts have [3.6] < 19.7 mag, suggesting that IRAC 
data to this depth would identify counterparts to ^ 60% of 
SPIRE sources. The remaining difference in the IRAC and 
WISE identification rates is likely to be due to the resolu- 
tions of the instruments - the WISE 3.4 /im beam is ~ 3 
times larger than IRAC - although Cosmic Variance may also 
contribute. We note that if instead the data were limited to 
[3.6] < 20.5, corresponding to the 5a point-source detection 
limit reached in a 120 sec integration with IRAC at 3.6 /im, 
the expected identification rate drops from 77% to ~ 70%. 

The results presented here are insensitive to the exact value 
of Qq. When using a Icr lower value for Qq (i.e. Qo = 0.84 at 
both 3.6 a nd 4 .5 /im), and comparing to the results presented 
in section [33] and table |2] the counterparts to 116 sources 
(95%) are unchanged. In this case there are four (3%) previ- 
ously identified sources that no longer have reliable counter- 
parts, and three (2%) previously unidentified SPIRE sources 
that now have R > 0.8 counterparts. 

3.4. Counterparts identified in IRAC color-magnitude space 

While we reliably identify 77% of the SPIRE sources with 
the LR method, the values of Qq indicate that ^ 90% of 
SPIRE counterparts are detected in IRAC. Figure |4] shows 
the IRAC color-magnitude diagram for galaxies in the H- 
ATLAS SDP field, with SPIRE counterparts identified with 
the LR method highlighted. These galaxies occupy a distinct 
region of IRAC color-magnitude space - they typically have 
[3.6] — [4.5] > — 0.4 and are brighter than the background 
population. We use this property to identify counterparts to 
the SPIRE sources that have both 3.6 and 4.5- /im IRAC cov- 
erage but no i? > 0.8 counterparts from the LR method. A 
similar analysis was performed by Biggs et al. (2011) at 3.6 
and 5.8-/im to identify to LABOCA 870- pm sources. 

We begin by determining the region of color-magnitude 
space in which there is a minimal chance of contamination by 
background sources (the ID region). A 20% contamination 
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TABLE 1 

Summary of IRAC counterparts to SPIRE sources 



Method 


A^spire" 




ID rate'= 


A''lRAc'' 


]\T e 
'multiple 


Contamination rate 


LR 


159 


123 


77% 


123 





1.9% 


Color-magnitude 


19 


13 


68% 


23 


7 


12.6% 


Total 


159 


136 


86% 


146 


7 


3.6% 



Note. — " Number of SPIRE sources considered. Number of SPIRE sources with at least one counterpart. Fraction of SPIRE sources with at least one 
identified counterpart. Number of IRAC counterparts identified. " Number of SPIRE sources with multiple IRAC counteiparts. 



limit is used because this equivalent to i? > 0.8 for the LR 
method. We demarcate the ID region with the simplest rea- 
sonable function - a single diagonal line. The placement of 
the line is determined by calculating the gradients and inter- 
cepts that would yield 20% contamination from background 
sources within 7.2" radius of the LR counterparts. The num- 
ber of LR counterparts returned is maximized to yield a best- 
fit gradient of 0.515 for a line that intercepts [3.6] = 16.0 mag 
at [3.6] — [4.5] = —2.142. This division is shown as a solid 
line in figure |4l IRAC galaxies that lie above this line and 
are within 7.2" of SPIRE sources that are otherwise uniden- 
tified, have > 80% probability of being physically associated 
with the SPIRE source, i.e., < 20% chance of being an unre- 
lated foreground or background galaxies, and are considered 
counterparts. 

There are 36 SPIRE sources that do not have LR coun- 
terparts. Of these 19 have both 3.6 and 4.5 //m data, and 
23 IRAC counterparts to 13 of these sources are identified 
with the color-magnitude method (seven have multiple coun- 
terparts). These sources are presented in table |2] The con- 
tamination rate for this method is 12.6% (approximately two 
IRAC galaxies), calculated by summing the probabilities of 
finding a background galaxy within 7.2" of a SPIRE source 
occupying a distinct region of color-magnitude space. 

In figure |4] we show the S'25o/'S'5oo color-flux plot for 
H-ATLAS SDP sources and consider whether sources with 
identified IRAC counterparts are representative of the whole 
SPIRE population. A 2D Kolmogorov-Smirnov (KS) test be- 
tween the background H-ATLAS SPIRE population and the 
sources with identified counterparts has p = 0. 172, suggesting 
that the two samples are drawn from similar, but not necessar- 
ily identical, parent populations. Considering the S'25o/'5'5oo 
color and 250 /im fluxes separately yields p — 0.710 and 
p = 0.289, respectively. Thus we conclude that the sources 
with identified IRAC counterparts are very similar to the 
whole SPIRE population, but may have a slightly different 
distribution of 250 /im flux. 

3.5. Astrometric offset 

To check whether there is any noticeable astrometric shift 
among the input catalogs for LR analysis, in figure|5]we plot 
the positional differences between the SPIRE sources and the 
corresponding IRAC counterparts. The median separation be- 
tween SPIRE and IRAC sources is 0.10". We also inves- 
tigate the overall astrometric reference frame difference be- 
tween IRAC and SDSS catalogs, using ~ 40% of the IRAC 
sample that is also identified with SDSS galaxies within in 2". 
The positional difference between IRAC and SDSS is 0.19 ± 
0. 17"and we found no systematic offset either in RA or Dec. 

Our IRAC data are complemented with VIKING Z, Y, J, H 
and Xg-band photometry, which is measured in 2" diam- 
eter apertures (Sutherland et al. in prep). We cross-match 
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Fig. 5. — The distribufion of ARA and ADec offsets between SPIRE 
sources and the con'esponding IRAC counterparts identified in this work. The 
astrometric median separation between SPIRE sources and IRAC is 0.10" 
which is consistent with the overall IRAC astrometric uncertainty tied to the 
SDSS astrometry. The three concentric circles have radii of 1, 2 and 3(jpos 
(7.2") from the SPIRE centroid. The overall offsets in RA and Dec show that 
there is no systematic offset between input SPIRE and IRAC sources. 

the IRAC (FWHM - 1.6") to the nearest neighbor VIKING 
(FWHM ^ 0.9") source within 2" search radius. Indeed, the 
mean offset between the IRAC and VIKING positions is 0.09 
± 0.50" and 0.28 ± 0.28" in RA and Dec, respectively. We 
find no statistically significant astrometric error between two 
catalogs. 

4. DISCUSSION 

As described in sections 13.31 and 13.41 we identified IRAC 
counterparts to 136 (86%) of the SPIRE sources. 123 of 
these are identified with the LR method and 23 counterparts 
to 13 SPIRE sources are identified with the color-magnitude 
method (Table[T]i. The contamination rate of 1.9% for the LR 
method is smaller than that from the color-magnitude method 
(12.6%) but corresponds to a similar number of IRAC galax- 
ies - two to three. Combining the results from the two iden- 
tification methods, we expect that the total false identification 
rate of the catalog presented in Table|2]is 3.6%, such that five 
to six of the IRAC counterparts presented in table|2]are false. 

4. 1 . Statistics of Identified Sources 

We next compare the counterparts identified in the IRAC 
data with results based on SDSS (Smith et al. 2011). Of the 
146 IRAC counterparts, 102 (70%) are undetected in SDSS, 
52 (36%) agree with the SDSS identification, and two (1%) 
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H -ATLAS 2244 (SPSS) H-ATLAS 2244 (3.6;^m) H-ATLAS 2244 (4.5)jm) 
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Fig. 6. — H-ATLAS ID-2244: an example of source that is identified using the LR method, but where the IRAC counteipart (R = 0.99 at both 3.6 and 4.5 ^m) 
is different from the SDSS counterpart (SDR2244, R = 0.97; Smith et al. 2011). From left-to-right we show a three-color SDSS image (g, r, and i bands), the 
3.6 fim, and the 4.5 fim data. The galaxies identified as the SDSS (SDP.2244) and IRAC (H-ATLAS ID-2244) counterparts to the SPIRE source are marked with 
an 'X' and a small circle, respectively. The large circle has 7.2" radius and encompasses the SPIRE 3o-pos area in which counterparts are identified. Contours 
show the SPIRE 250 /im emission at 5, 7, 9, Hit levels. In the absence of high-resolution sub-mm imaging we cannot determine whether this SPIRE source is 
the results of blended emission from the two identified galaxies, or whether one of those counterparts is a chance association. 



.^Tt^-^— iiM^Ii^^-^-^— I ■- ■ — H-ATLAS 891 C3.6/im) H-ATLAS 891 {4.5/im) 
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Fig. 7. — H-ATLAS ID-891: an example of counterpart identification using the color-magnitude method. Symbols are as figure|6] There is no LR counterpail 
in the SDSS r-band (Smith et al. 201 1) or with IRAC at 3.6 or 4.5 ^m. However, there are two galaxies that are both identified as SPIRE counterparts on the 
basis of their IRAC colors and magnitudes. These two galaxies have R. = 0.41 and 0.47 at 3.6 /im, and R = 0.44 and 0.50 at 4.5 fj,m and neither is detected 
in SDSS. It is probable that the SPIRE source is a blend of the emission from these two IRAC galaxies although, in the absence of spectroscopic redshifts it is 
unclear whether the configuration is caused by a line-of-sight alignment or an interaction between the two IRAC galaxies. 



H-ATLAS 827 (SDSS) H-ATLAS 827 (3.6/im) H-ATLAS 827 (4.5/zm) 
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Fig. 8. — H-ATLAS ID-827 has four counteiparts identified with the IRAC color-magnitude method and is the most complex system in our sample. The 
westernmost and easternmost galaxies are detected by SDSS (but are not identified as counterparts in these data; Smith et al. 2011), and have photometric 
redshifts of z = 0.07 ± 0.01 and z = 0.20 ± 0.03, respectively. The two fainter counterparts are not detected in SDSS but have [3.6] — [4.5] > 0, indicating that 
they may be at z > 1.4 (Papovich 2008). We conclude that the SPIRE source is most likely to be compiised of blended emission from the four galaxies. Symbols 
are as figure |6] 
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Fig. 9.— The distribution of IRAC [3.6] - [4.5] color for SPIRE coun- 
terparts witli pliotometry at botli wavelengtlis. For tlie 47 sources with 
SDSS spectroscopic or photometric redshifts we also plot the [3.6] — [4.5] 
color against redshift. The 57 SPIRE-IRAC identifications without reli- 
able redshifts are either faint or undetected in SDSS and 43 of them have 
[3.6] — [4.5] > 0, an indicator of high redshift source (z > 1.4; Papovich 
2008). We expect that more than half of the sub-mm bright sources reside at 
2 > 1.4. For reference the color-redshift tracks of Arp 220 (Silva et al. 1998), 
Mrk23 1 (Berta 2006) and the average sub-mm galaxy SED (Wardlow et al. 
201 1) are shown. 



are different galaxies to the SDSS counterparts. On the basis 
of the false-identification rates of the two studies (4.2% for 
SDSS and 3.6% for IRAC) approximately four counterparts 
are expected to disagree between the two surveys, which is 
consistent with that observed here. 

The two sources with different counterparts identified in 
IRAC and SDSS are H-ATLAS ID-2244 and 6962. We show 
one of these sources, H-ATLAS ID-2244, as an example in 
figure |6] In this case the SDSS (SDR2244) and IRAC (H- 
ATLAS ID-2244) counterparts are separated by 3.2" and both 
have a < 3% probability of being false identifications: R = 
0.97 in SDSS and R = 0.99 at both 3.6 and 4.5 /im. In the 
absence of additional data, such as sub-mm or radio interfer- 
ometry, the true nature of this source is unclear It is possible 
that one of the counterparts is a chance association, although 
we cannot say which one. It is also possible that the SPIRE 
source is comprised of a blend of emission from both coun- 
terparts in a single SPIRE beam. 

Blending is also an important consideration for counterparts 
that are identified with the color-magnitude method. Seven 
of the 13 SPIRE sources identified with this method have 
multiple IRAC counterparts; conversely the LR method has 
none. This apparent disparity is not surprising because the 
LR method implicitly assumes that each SPIRE source has a 
single IRAC counterpart. Where there are N multiple coun- 
terparts contributing equally to the sub-mm flux the average 
R cannot exceed 1/iV, and will typically be ^ 1/iV. Due to 
the large SPIRE beam (FWHM = 1 8. 1 " at 250 ^m) it is likely 
than in at least a few cases a single SPIRE source may be 
composed of blended emission from multiple galaxies or from 
galaxy interactions. The color-magnitude method does not 
consider the presence of other IRAC sources and is therefore 
not biased against multiple counterparts. In addition, beyond 
requiring that counterparts to lie within the SPIRE counterpart 
search radius, the color-magnitude method does not consider 
the separation between the SPIRE and IRAC centroids and 
therefore it is not biased against wide-separation counterparts. 



However, it does require both 3.6 and 4.5 /zm fluxes and it as- 
sumes that all SPIRE sources have a similar color-magnitude 
distribution at 3.6 and 4.5 /im. 

The seven SPIRE sources in the sample that have multiple 
IRAC counterparts are likely to be a combination of blended 
SPIRE sources and individual sub-mm galaxies associated 
with multiple interacting IRAC sources. H-ATLAS ID-891 is 
shown as an example in figure|2l In this case all the detected 
SDSS galaxies are > 3.5" away from the SPIRE centroid and 
all have R < 0.8 (Smith et al. 201 1). There are two potential 
counterparts in the IRAC data that are not detected in SDSS. 
Using the LR method alone these two nearby sources have re- 
liability of i? = 0.41 and 0.44 at 3.6 and R = 0.47 and 
0.50 at 4.5 /im. However, both galaxies have [3.6] and [4.5] 
that place them in the region of color-magnitude space with < 
20% probability of being random associations with the SPIRE 
emission and as such they are both identified as counterparts. 
The two galaxies have similar colors, with ([3.6] — [4.5]) — 
0.38 and 0.43 mag, consistent with a redshift of z > 1.4 (Pa- 
povich 2008). The two IRAC sources are separated by 2.9" 
(24.9 kpc at 2 =2, 23.0 kpc at z =3. Furthermore, the 
SPIRE source is bright (55oo = 84 ± 9mJy) and red in the 
sub-mm bands, with a rising spectrum from 250 to 500 i-im, 
indicating a sub-mm photometric redshift of z > 3 (e.g. Lapi 
et al. 201 1), although we caution that the dust temperature and 
redshift are degenerate in the sub-mm bands. It is likely that 
H-ATLAS ID-891 is a blend of emission from the two identi- 
fied counterparts and it is possible that the sub-mm emission 
is the result of a merger or interaction of two gas-rich galaxies 
(e.g. Aravena et al. 2010b). Spectroscopic data are required 
to verify this scenario, and we identify ID-891 as a target for 
additional follow-up, especially to identify the exact nature of 
the sub-mm emission. 

H-ATLAS ID-827, shown in figure [8] is the most complex 
system in our sample. In this case there are four IRAC coun- 
terparts identified on the basis of their colors and magnitudes. 
The four galaxies have a wide range of [3.6] with values from 
17.8 mag to 20.2 mag. The counterparts do not appeared to 
be clustered, as may be expected for gravitational lensing or 
a multi-component interaction, and each is separated from its 
nearest neighbor by ^ 4". The two brightest galaxies are de- 
tected by SDSS, but were not identified as SPIRE counter- 
parts by Smith et al. (2011); they have photometric redshifts 
of z = 0.07 ± 0.01 and z = 0.20 ± 0.03 for the western- 
most and easternmost, respectively. The two fainter counter- 
parts are not detected in SDSS and both have [3.6] — [4.5] > 
0, indicating that they are likely to be at z > 1.4 (Papovich 
2008). The morphology, astrometry and available redshift in- 
formation suggests that gravitational lensing in unlikely and 
that H-ATLAS ID-827 is most likely to be a blended source, 
which is comprised of two low-redshift and two high-redshift 
components. 

4.2. Redshifts of Identified Sources 

Knowledge of the redshift distribution of SPIRE sources 
is central to understanding their role in the Universe. How- 
ever, only 50% of the IRAC counterparts are detected in SDSS 
and 74% have VIKING photometry. Therefore, any redshift 
distribution derived from the optical data will be biased to 
low redshifts where the detection rates in these surveys are 
high. Instead, in figure|9]we plot the [3.6] — [4.5] color dis- 
tribution, and [3.6] — [4.5] against redshift for the counter- 
parts with spectroscopic or photometric redshift from SDSS. 
This plot includes the 104 IRAC counterparts that have both 
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Fig. 10.— IRAC and VIKING color-color diagrams. Top: [Ks], [3.6] 
and [4.5]. Bottom: [z'], [Ks] and [3.6]. We highlight SPIRE counterparts 
and background sources that are within 7.2" of SPIRE centroids. Tracks for 
Aip 220 (Silva et al. 1998), Mrk 231 (Berta 2006) and the average 870- fim 
selected sub-mm galaxies (Wardlow et al. 2011) are shown, and an Ay = 
1 mag reddening vector and the median error bai' are plotted in the lower left- 
hand comer. The SPIRE counterparts ai'e slightly bluer than both these tracks 
suggesting that they may have less dust reddening and Ay values that are 0.5 
to 1 magnitudes smaller than the SED templates. 

3.6 and 4.5 i-im coverage, 47 of which have an optical red- 
shift. All of the optical redshift are z < 0.8; 32 (68%) of 
these galaxies have ([3.6] - [4.5]) < 0, while 15 (32%) have 
([3.6] — [4.5]) > 0. The 57 identifications without reliable red- 
shifts from the optical are undetected in SDSS and 43 of them 
have ([3.6] — [4.5]) > 0, a crude indicator that z > 1.4 (e.g. 
Papovich 2008). Thus, up to - 40% of the identified SPIRE 
population may lie at high redshifts. Deep spectroscopy and 
photometry, particularly at near-lR wavelengths, is required 
to determine more precise redshift information. 

We also investigate the redshift distribution of the SPIRE 
sources by calculating sub-mm photometric redshifts from the 
fluxes of the three SPIRE bands. We employ a template fitting 
method which uses minimization, comparing the observed 
fluxes to the SEDs of Arp 220 (Silva et al. 1998), Mrk 231 
(Berta 2006). We caution that the sub-mm photometric red- 
shifts are limited by the absence of longer wavelength data, 
and the degeneracy between redshift and dust temperature. 
Furthermore, flux boosting is an important consideration and 
will add significant errors to the fluxes of the faintest sources. 
Due to these uncertainties, we do not consider the sub-mm 



Fig . 11 . — Top: Dust temperature as a function of the infrai'ed luminosity 
LiR (8 - 1100 ^tm) for H-ATLAS/IRAC galaxies with SDSS redshifts. Lir is 
obtained by fitting the SPIRE flux with /3 = 1.5 in isothermal modified black- 
body model (equation |9). EiTor bai's ai'e median la standard deviation from 
the best-fit models with dust temperature (Tj) and total IR luminosities Lir 
as free parameters. Bottom: Dust temperature as a function of the redshift 
for H-ATLAS/IRAC galaxies. See text for the comparison galaxy samples 
plotted here. 

photometric redshifts on a source-by-source basis, but instead 
examine the overall distribution. 

This sub-mm photometric redshift distribution peaks at z '-^ 
2 with a tail out to 2; 4, and suggests that the majority of the 
IRAC counterparts without SDSS redshifts lie at z ^ 2. This 
result is consistent with the photometric redshifts of the whole 
SPIRE population (e.g. Amblard et al. 2010; Lapi et al. 201 1) 
and the redshift distribution of 850 i-im sources (e.g. Chapman 
et al. 2005; Wardlow et al. 201 1). Now that the counterparts to 
a significant majority of SPIRE sources are known, follow-up 
spectroscopic campaigns to establish the redshifts of all the 
sources are feasible. 

4.3. Properties of Identified Sources 

Near-IR Colors: In figure [lO] we show IRAC and VIKING 
color-color diagrams, highlighting identified SPIRE counter- 
parts in comparison with background sources that are within 
7.2" of the SPIRE centroids. As we show in figure 2] the 
SPIRE counterparts have distinct [3.6] — [4.5] color compared 
to the background sources, but the two populations are indis- 
tinguishable in [/Vs] — [3.6] and [z] — [Kj^. We also com- 
pare to the redshift tracks of Arp 220 (Silva et al. 1998), 
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Fig. 12. — Top: Specific star formation rate (sSFR) for our IRAC- 
identified SPIRE galaxies as a function of the redshift. Filled circles and 
triangles represent H-ATLAS galaxies in stellai' mass bins lO'-* **- lO^*^'* 
Mq and lO^'' *- 10^^'^ Mq, respectively. Our data are compared to Perez- 
Gozalez et al. (2008) and Damen et al. (2009). Upper (lower) two lines are 
for the smaller (larger) mass bin. The eiTor bar on the top left shows the 
average 68% confidence range on iog(sSFR). Bottom: Star formation rate 
as a function of the stellar mass. We plot results from a previous study for 
comparison: PACS-COSMOS and PACS-GOODS South (Rodighiero et al. 
2011). We also show the mean SFR and stellar mass correlation at z = 2 
(Daddi et al. 2007), 1 and (Elbaz et al. 2007). The majority of the H- 
ATLAS/IRAC galaxies fall between 2 = ~ 2 in SFR-A/* relation showing 
a slight excess of star formation rate compared with populations at similar' 
redshifts, but selected at optical and near-IR wavelengths. 



Mrk 23 1 (Berta 2006) and the average 870- fim selected sub- 
mm galaxy (Wardlow et al.2011). The SPIRE counterparts 
are slightly bluer than both these tracks suggesting that they 
may have less dust reddening and Ay values that are 0.5 to 1 
magnitudes smaller than the SED templates. 

IR Luminosities and Dust Temperatures: To understand 
properties of the identified H-ATLAS/IRAC galaxies associ- 
ated with SDSS redshifts, we first consider a SED analysis of 
the SPIRE data making use of the redshifts based on the IRAC 
identification. First we consider the far-lR/sub-mm portion of 
the SED with 3-band SPIRE fluxes at 250, 350, and 500 /im 
and reproduce the analysis of Amblard et al. (2010) where the 
dust temperature (Td) and IR luminosities Lir were studied. 
We make use of a modified black-body with isothermal dust 
temperature to describe the dust emission, in which flux can 
be written as 

^ocM'+^e('^^/'^^--i). (9) 



To be consistent with previous other estimates from the litera- 
ture we fix /3 = 1.5 and fit to the SPIRE fluxes to establish the 
dust temperature and the overall normalization to the SED. 
The IR luminosity is estimated from the best-fit SED over the 
wavelength range of 8 to 1 100 /im. Since we only fit to the 
dust emission at wavelengths greater than 250 fim our lumi- 
nosities are likely low by a factor of 2 to 3 if there are any 
warm dust, heated by AGNs, present in these galaxies. How- 
ever, our overall uncertainty on the IR luminosity estimates is 
at least a factor of 5. For the sample as a whole we find an 
average dust temperature of 26 ± 9 K, which compares well 
with the average dust temperature of 27 ± 8 K for z < O.l 
SPIRE galaxies based on SDSS identifications only (Amblard 
et al. 2010). 

We summarize our results related to the FIR SED analysis 
in figure [TT| where we show the Lir - relation and the Lir 
as a function of the redshift, determined from the reliable H- 
ATLAS IDs with either SDSS spectroscopic or photometric 
redshift. We also compare the H-ATLAS dust temperatures 
and IR luminosities with sub-mm bright galaxy samples in 
the literature. Samples include the sources in BLAST with 
COMBO- 17 (Wolf et al. 2004) or a SWIRE photometric red- 
shift (Rowan-Robinson et al. 2008) and Spitzer-MlPS 70 and 
160 /im fluxes (Dye et al. 2009), local ULIRGS observed with 
SCUBA at 450 and 850 /tm (Clements et al. 2010), and local 
IRAS-selected galaxies with 60 and 100 /im fluxes along with 
SCUBA 850 /im (Dunne et al. 2000). 

In figure [Til we find that the SPIRE-selected galaxies in H- 
ATLAS with IRAC -based identifications for redshifts span 
the luminosity range of 10^" to 10^^ Lq from sub-LlRG 
luminosities to ULIRG conditions. These luminosities and 
the ranges are comparable to BLAST-detected sources, while 
they span lower than the local IRAS-selected galaxies, which 
tend to be at lower redshifts z < 0.05 and have higher lumi- 
nosities; such a difference is understandable since IRAS is an 
all-sky shallow survey and detects primarily the rare, bright 
galaxies in the near-by universe. 

Stellar Masses and Star-Formation Rates: With the cross- 
identification we can expand the SED analysis over 3 orders of 
magnitude in wavelength from optical to sub-mm. For optical 
and near-lR data we make use of SDSS, UKIDSS (Lawrence 
et al. 2007), VIKING (Sutherland et al. in prep), and IRAC 
3.6 and 4.5 /im fluxes. The analysis of full optical to sub- 
mm SEDs is similar in spirit to Dunne et al. (2011) where 
the SEDs of SPIRE sources with reliable identifications using 
SDSS data were analyzed. We make use of the MAGPHYS 
SED modelling code (da Cunha, Chai'lot & Elbaz 2008) for 
this work. The SED modelling in the code is done in a energy- 
balance manner such that the absorbed light in shorter UV 
and optical wavelengths is accounted for by the thermal re- 
radiation at far-infrared and sub-mm wavelengths. For the H- 
ATLAS/IRAC galaxy sample we have 5 optical fluxes from 
SDSS (it, g, r, i, z), 5 near-IR fluxes from UKIDSSh-VIKING 
(Z, Y, J, H, Ks), and 2 IRAC channels (3.6, 4.5 /im), in addi- 
tion to 3 SPIRE bands. Redshift information of the sources 
are from SDSS spectroscopic observations (z =0 - 0.8). 

For comparative work with other sub-mm galaxy samples 
we focus on results related to stellar mass Mi, and star- 
formation rate (SFR), as given by MAGPHYS best-fit mod- 
els. In figure[T2](top figure) we plot the specific star-formation 
rate (sSFR) given as SFR/A/^, vs. redshift with our sample 
sub-divided to two stellar mass bins. We find that the sSFR 
of IRAC IDs marginally show an anticorrelation with galaxy 
stellar mass. Such a behavior may be related to the downsiz- 
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ing scenario (Cowie et al. 1996). 

However, when compared to Damen et al. (2009) and 
Perez-Gozalez et al. (2008), IRAC IDs shows 0.4 - 0.5 dex 
higher sSFRs. The Damen et al. (2009) galaxy sample 
is an IRAC-selected sample of galaxies out to z ~ 1.8 in 
ECDFS with photometry supplemented with optical, near-IR, 
and MIPS 24 fim fluxes. The Perez-Gozalez et al. (2008) sam- 
ple is selected with IRAC at 3.6 and 4.5 i-im in a variety of 
fields. While our SPlRE-selected sample has IRAC identifi- 
cations, they are not likely to be the majority of the galax- 
ies in both these studies. This is also clear from figure |4] 
where we show that most of the identified SPIRE sources are 
distributed differently from the background population in the 
color-magnitude space. They occupy the top-end of the star- 
formation in galaxies selected at near-lR wavelengths. 

In Fig. 12 (bottom) we plot the SFR vs. the stellar mass 
of our sample. The SFR of normal local star-forming galax- 
ies at z ~ 1 is known to correlate strongly with the stellar 
mass and form the so-called star- formation main-sequence 
(Brinchmann et al. 2004; Salim et al. 2007; Peng et al. 2010). 
Such a correlation is also observed at higher redshifts, albeit 
with a different normalization (e.g. Daddi et al. 2007; Elbaz 
et al. 2007; Noeske et al. 2007; Pannella et al. 2009; Daddi et 
al. 2009; Gonzalez et al. 2010; Rodighiero et al. 2010; Karim 
et al. 201 1). However, the typical dispersion of correlation is 
known to be ^ 0.3 dex over a wide range of redshifts z = 
~2. 

Local (U)LlRGs and high-redshift SMGs have SFRs in 
excess of the main-sequence and are defined as starburst 
galaxies (e.g. Elbaz et al. 2007, 2011). These galaxies 
generally occupy a region that is ^ xlO above the main- 
sequence. However, there is also evidence that some high- 
redshift (U)LlRGs and SMGs do not have enhanced star- 
formation rates relative to their stellar mass and that the star- 
burst fraction may decrease at high redshifts (e.g. Tacconi et 
al. 2008; Rodighiero et al. 2011). This suggests changes to 
the mode of star-formation in (U)LlRGs between z = and 
z = 2 away from the starburst mechanism (e.g. Daddi et al. 
2010b; Genzel et al. 2010; Elbaz et al. 2011; Krumholz et 
al. 2012; Melbourne et al. 2012). In our work, - 50% of 
the subsample of IRAC-identified SPIRE galaxies (z < 1) lie 
near the starburst region. Thus, it is possible that half of the 
IRAC-identified SPIRE galaxies with SDSS redshifts ai-e in 
fact starbursts that lie off the main-sequence. 

5. SUMMARY 

We have identified Spitzer-IKAC counterparts to sources 
selected with //erac/ieZ-SPlRE at 250 /im in the H- ATLAS 
survey. Among 159 SPIRE centroids, we found 123 reliable 
IRAC counterparts using a likelihood ratio analysis. The iden- 
tified SPIRE sources are distributed differently in IRAC color- 
magnitude space compared to the field population. There- 
fore, we made an additional selection based on the IRAC 
ID locality to yield a further 23 counterparts to 13 SPIRE 
sources. Seven SPIRE sources have multiple IRAC counter- 
parts. These are likely to be due to blended emission in the 
SPIRE beam. 

In total 146 reliable IRAC counterparts to 136 SPIRE 
sources were identified, including seven that have more than 
one IRAC counterpart. The identification rate of 86% is 



higher than that of wide-field ground-based optical and near- 
IR imaging of Herschel fields. The galaxies with unknown 
redshifts and that are not detected in SDSS and VIKING 
imaging data have SPIRE colors indicative of high redshift 
sources and [3.6] — [4.5] > 0, suggesting that they are likely 
to be at z > 1.4. We estimate the 40% of SPIRE sources 
lie at high redshifts, although the exact redshift distribution of 
SPIRE sources remains elusive. The counterparts presented 
here can now be pursued for followup data to further inves- 
tigate the nature of SPIRE galaxies. The majority of our de- 
tected galaxies with sub-LIRG to LIRG luminosities are not 
intense starbursting galaxies in the local universe, though they 
have above average specific star-formation rates. 



Herschel-ATLAS is a project which uses Herschel, an 
ESA space observatory with science instruments provided by 
European-led Principal Investigator consortia and with im- 
portant participation from NASA. The H-ATLAS website is 
http://www.h-atlas.org/. This work is based in part on ob- 
servations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Insti- 
tute of Technology under a contract with NASA. Support for 
this work was provided by NASA through an award issued by 
JPL/Caltech. US participants in H-ATLAS also acknowledge 
support from the NASA Herschel Science Center through a 
contract from JPL/Caltech. 

Funding for the SDSS and SDSS-II has been provided by 
the Alfred P. Sloan Foundation, the Participating Institutions, 
the National Science Foundation, the U.S. Department of En- 
ergy, the National Aeronautics and Space Administration, the 
Japanese Monbukagakusho, the Max Planck Society, and the 
Higher Education Funding Council for England. The SDSS 
Web Site is http://www.sdss.org/. 

The SDSS is managed by the Astrophysical Research Con- 
sortium for the Participating Institutions. The Participating 
Institutions are the American Museum of Natural History, As- 
trophysical Institute Potsdam, University of Basel, Univer- 
sity of Cambridge, Case Western Reserve University, Uni- 
versity of Chicago, Drexel University, Fermilab, the Institute 
for Advanced Study, the Japan Participation Group, Johns 
Hopkins University, the Joint Institute for Nuclear Astro- 
physics, the Kavli Institute for Particle Astrophysics and Cos- 
mology, the Korean Scientist Group, the Chinese Academy 
of Sciences (LAMOST), Los Alamos National Laboratory, 
the Max-Planck-Institute for Astronomy (MPIA), the Max- 
Planck-Institute for Astrophysics (MPA), New Mexico State 
University, Ohio State University, University of Pittsburgh, 
University of Portsmouth, Princeton University, the United 
States Naval Observatory, and the University of Washington. 

This work used data from VISTA at the ESO Paranal Ob- 
servatory, programme 179.A-2004. VISTA is an ESO near- 
infrared survey telescope in Chile, conceived and developed 
by a consortium of universities in the United Kingdom, led by 
Queen Mary University, London. 

We thank Giulia Rodighiero for providing us electronic 
tables of COSMOS and GOODS SED fitting results. The 
UCI group acknowledges support from NSF CAREER AST- 
0645427, NASA NNX10AD42G, and an award from Cal- 
tech/JPL for US participation in Herschel-ATLAS. 



REFERENCES 



Abazajian, K. N. et al. 2009, ApJS, 182, 543 



Amblard, A. et al. 2010, A&A, 518, L9 



13 



Aiavena, M. et al. 2010a, ApJ, 708, L36 
Aiavena, M. et al. 2010b, ApJ, 719, 15 
Ashby, M. L. N. et al. 2009, ApJ, 701, 428 
Austermann, J. E. et al. 2010, MNRAS, 401, 160 
Bai-mby, P. et al. 2008, ApJS, 177, 431 
Berta S. 2006, PASP, 118, 754 
Bertin, E. & Arnouts, S. 1996, A&AS, 1 17, 393 
Biggs, A. D. et al. 2011, MNRAS, 413, 2314 

Blain, A. W., Small, 1., Ivison, R. J., Kneib, J.-P & Prayer, D. T. 2002, 

Physics Reports, 369, 111 
Bond. N., et al. 2012, ApJL, 750, 18L 
Borys, C, et al. 2004, MNRAS, 355, 485 
Brandon, K., et al. 2006, ApJ, 641. 140 
Brinchmann, J., et al. 2004. MNRAS, 351, 1151 
Brusa, M., et al. 2007, ApJS, 172, 353 
Budavari, T. & Szalay. A. 2008. ApJ, 679, 301 
Chapin, E.. et al. 2009, MNRAS. 398. 1793 

Chapman, S. C, Blain, A. W., Ivison, R. J., & Small, I. 2003, Nature, 422, 
695 

Chapman, S. C, Blain, A. W., Small, I., & Ivison, R. J. 2005, ApJ, 622, 772 

Chen, C.-C, et al. 2011, ApJ, 733, 64 

Ciliegi, R, et al. 2003, A&A, 398, 901 

Clements, D. L., et al. 2008, MNRAS, 387, 247 

Clements, D. L., Dunne, L., & Bales, S. 2010, MNRAS, 403, 274 

Condon, J J., et al. 1998, AJ, 1 15, 1693 

Coppin, K., et al. 2006, MNRAS, 372, 1621 

Cowie, L. L., et al. 1996, AJ. 112. 839 

Cowie, L. L., et al. 2009, ApJ, 697, L122 

da Cunha, E., Chariot, S., & Elbaz, D. 2008, MNRAS, 388, 1595 

Daddi, E., et al. 2007, ApJ, 670, 156 

Daddi, E., et al. 2009, ApJ, 694, 1517 

Daddi, E., et al. 2010, ApJ, 714, 1 18 

Damen, M., et al. 2009, ApJ, 690, 937 

Downes, D., et al. 1999, A&A, 347, 809 

Driver, S., et al. 2011, MNRAS, 413, 971 

Dunlop, J. S., et al. 2010, MNRAS, 408, 2022 

Dunne, L., Eales, S., Edmunds, M., et al. 2000, MNRAS, 315, 115 

Dunne, L., et al. 201 1, MNRAS, 417, 1510 

Dwek, E., et al. 1998, ApJ, 508, 106 

Dye, S., et al. 2008, MNRAS, 635, 1 107 

Dye, S., et al. 2009, ApJ, 703, 285 

Eales, S., etal. 1999, ApJ, 515, 518 

Eales, S., et al. 2010, PASP 122, 499 

Eisenhardt, R R., et al. 2004, ApJS, 154, 48 

Elbaz, D., et al. 2007, A&A, 468, 33 

Elbaz, D., et al. 2011, A&A, 533, 119 

Fazio, G. G., et al. 2004a, ApJS, 154, 10 

Fazio, G. G., et al. 2004b, ApJS, 154, 39 

Fixsen, D. J., Dwek, E., Mather, J. C, Bennett, C. L., & Shafer, R. A. 1998, 
ApJ, 508, 123 

Fleuren, S., et al. 2012, MNRAS, accepted (arXiv:1202.3891) 

Garrett, M. A. 2002. A&A, 384, 19 

Gear, W. K., et al. 2000, MNRAS, 316, 51 

Genzel, R., et al. 2010, MNRAS, 407. 209 

Gonzalez, V., et al. 2010, 713, 115 

Gorjian, V., et al. 2008, ApJ, 679, 1040 

Griffin, M. J., et al. 2010, A&A, 518, L3 

Hatsukade, B., et al. 2010, ApJ, 71 1, 974 



Hopwood, R., et al. 201 1, ApJ, 728, L4 

Hughes, D. H., et al. 1998, Nature, 394, 241 

Ibar, E., et al. 2010, MNRAS, 409, 38 

lono, D., et al. 2006, ApJL, 640, LI 

Ivison, R. J., et al. 1998, ApJ, 494, 21 1 

Ivison, R. J., et al. 2002, MNRAS, 337, 1 

Ivison, R. J., et al. 2004, ApJS, 154, 124 

Ivison, R. J., et al. 2007, MNRAS, 380, 199 

Ivison, R. J., et al. 2008, MNRAS, 390, 1117 

Karim, A., et al. 2011, ApJ, 730, 61 

Krumholz, M., et al. 2012, ApJ, 745, 69 

Lapi, A., etal. 2011,742, 24 

Lawrence, A., et al. 2007, MNRAS, 379, 1599 

Lutz, D., et al. 201 1, A&A, 532, A90 

Makovoz, D., & Marleau, F. R. 2005, PASP 1 17, 1 1 13 

Matsuda, Y., et al. 2007, ApJ, 667, 667 

Melbourne, J., et al. 2012, AJ, 143, 125 

Negrello, M., et al. 2010, Science, 330, 800 

Noeske, K. G., et al. 2007, ApJL, 660, L43 

Pannella, M., et al. 2009. ApJ, 701, 787 

Papovich, C. 2008, ApJ, 676, 206 

Pascale, E., et al. 2011, MNRAS, 415, 911 

Peng, Y.-J.. et al. 2010, ApJ, 721, 193 

Perez-Gonzalez, R G., et al. 2008. ApJ, 675, 234 

Pilbratt, J., et al. 2010, A&A, 518, LI 

Poglitsch, A., et al. 2010, A&A, 518, L2 

Pope, A., et al. 2006, MNRAS, 370, 1185 

Puget, J. L., et al. 1996. A&A, 308, L5 

Rigby, E. E., et al. 201 1, MNRAS, 415, 2336 

Rodighiero, G., et al. 2010, A&A, 518, 25 
Rodighiero, G., et al. 201 1, ApJL, 739L, 40 

Rowan-Robinson, M., Babbedge, T., Oliver, S., et al. 2008, MNRAS, 386, 
697 

Salim, S., et al. 2007, ApJS, 173, 267S 
Schneider, D. R, et al 2007, AJ, 134, 102 
Scott, S. E., et al. 2002, MNRAS, 331, 817 
Scoville, N., et al. 2007, ApJS, 172, 1 

Silva, L., Granato, G. L., Bressan, A., & Danese, L. 1998, ApJ, 509, 103 

Small, I., Ivison, R. J., Blain, A. W. 1997, ApJ, 490, L5 

Smolcic, v., et al. 2012a, ApJS, 200, 10 

Smolcic, v., et al. 2012b, A&A, accepted (arXiv:1205.6470) 

Smith, D. J. B., et al. 201 1, MNRAS, 416, 857 

Surace, J. A., Shupe, D. L., Fang, F., Lonsdale, C. J., & Gonzalez-Solares et 

al. 2005. http://swire.ipac.caltech.edu/swire/astronomers/pubhcations 
Sutherland, W. & Saunders, W. 1992, MNRAS, 259, 413 
Swinbank, A. M., Smail, I., Longmore, S., et al. 2010. Nature, 464, 733 
Tacconi, L. J., et al. 2008, ApJ, 680, 246 
Tamura, Y, et al. 2010, ApJ, 724, 1270 
Thompson, M., et al. 2010, A&A, 518, L134 
Wang, W.-H., et al. 2011, ApJ, 726, L18 
Wardlow, J. L., et al. 2011, MNRAS, 415, 1479 

Wolf, C, Meisenheimer, K., Kleinheinrich, M., et al. 2004, A&A, 421, 913 
Younger, J. D., et al. 2007, ApJ, 671, 1531 
Younger, J. D., et al. 2008, MNRAS, 387, 707 
Younger, J. D., et al. 2009, ApJ, 704, 803 
York, D. G., et al. 2000, AJ, 120, 1579 



APPENDIX 



TABLE 2 

Combined photometry for IRAC IDs. 
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01 


18 


02 


± 





01 


17 


75 


± 





01 


17 


94 


± 





01 












46 


36 


± 


6 


47 


16 


00 


± 


7 


20 


-2 


.30 


± 9 


.00 





998 









20 


± 





04 


4 


56 



Note. - * ID with color-magnitude method. " Coordinates are the position of IRAC counterparts. ^ Photometry from the VIKING survey {Z, Y, J, H, Kg) is extracted in 2" diameter apertures 
(Sutherland et al. in prep). The reported SPIRE 350 and 500 /im fluxes are measured at the postions of 250 /im sources; no SNR cut is applied. R is reliability of the counteiparts calculated separately at 3.6 and 
4.5 fim from the LR method (section lsTl . Redshifts are from the SDSS catalogue for sources detected in that survey, f The separation between the centroid of the SPIRE emission and the IRAC counterpart. 



TABLE 2 - continued. 



ID 



(J20QQ) 

"08^41^38 

■^08^49^79 
"08^50?16 
"08^03?57 
"08^55?59 
"09^30?22 
"08^03?82 
"09^05? 49 
"08^ 46? 01 
"08^45?88 
"09^43? 23 
"08^20?Q0 
"08^02?92 
"09h22?19 
"03^33? 53 
"08^00?51 
"08^39?11 
"07h46?3S 
"08^32?85 
"08^28?96 
"08^43^61 
"08^03?34 
"08^36?e2 
"09^40? 28 
"09^21?78 
"08^32?24 
"03^23?60 
"03^ 17? 75 
"09^01?77 
"08^58?82 
"08^01?44 
"09^03?34 
"09^22? 24 
"09^01?87 
"08^12?06 
"08^21?24 
"09^35?26 
"09^18^30 
"09^18?44 
"09^07?08 
"09^18?11 
"07^44? 70 
"09^02?05 
"09^02?39 
"09^01?94 
"09^04?37 
"08^51?ee 



(J2000) 



4.5 /^m 

(mag) 



•5250'' 
(mJy) 



'^350'^ 
(mJy) 



■5500'' 
(mJy) 



(3.6^m) (4.5^m) 



Sopiiralion/ 



J090925. 
J091257. 
J090841. 
J090752. 
J090850. 
J090850. 
J090803. 
J090855. 
J09093O. 
J090803. 
J090905. 
J090846. 
J090846. 
J090943. 
J090819. 
J091302. 
J090922. 
J090333. 
J090800. 
J090839. 
J090746. 
J090832. 
J090829. 
J090843. 
J091303. 
J090836. 
J090940. 
J090921 
J090832. 
J090323. 
J0913!7. 
J090901. 
J090858. 
J09100 
J090903. 
J090921 
J090901 
J090811 
J090821 
J090935. 
J090918. 
J090918. 
J091306. 
J091318. 
J090744. 
J090902. 
J090902. 
J090902. 
J091304. 
J09085 



2+003224 
8-005508 
3+002005 
3+002100 
.0+004309 
.0+004309 
8+002250 
5+002808 
2+002755 
.7+002921 
3+001525 
.0+004339 
.0+004339 
.0+004322 
.6+003259 
.7-004618 
.4+002715 
3+004746 
5+002457 
.4+004107 
.7+002148 
.7+002406 
1+001556 
.7+001437 
.6-004855 
.4+002948 
3+002939 
.9+002556 
2+001938 
.9+004620 
.7-004621 
5+003107 
8+003158 
.4+004024 
5+00203 1 
.9+004307 
.9+004217 
.9+003410 
.6+002700 
1+002224 
3+003409 
3+003409 
.9-004719 
1-005409 
.7+002005 
.0+001936 
.0+001936 
.0+001936 
1-005141 
.6+003823 



2281 

2343 

2432 

2437 

2459* 

2459* 

2549 

2565 

2680 

2715 

2773 

2793* 

2793* 

2796 

2866 

2986 

3043 

3056 

3084 

3113 

3161 

3205 

3242 

3251 

3285 

3304 

3318 

3457 

3573 

3574 

3583 

3616 

3714 

3792 

3825 

3974 

4113 

4185 

4344 

4352 

4366* 

4366* 

4404 

4520 

4524 

4587* 

4587* 

4587* 

4658 

4693 



+ 00°32''26' 
-00°55'02' 
+00°20'06' 
+00°2l'00' 
+00°43'07' 
+00°43'09' 
+00°22''51 

+oo°28'oe' 

+00°27'55' 
+00°29'22' 
+00° 15^23' 
+00°43'36' 
+00°43'39' 
+00°43''22' 
+00° 33^02' 
-00°46'l9' 
+00° 27' 15' 
+00°47'48' 
+00°24'57' 
+00°4l'06' 
+00°21''47' 
+00°24'07' 
+00°15'5S' 
+00° 14^38' 
-00°48'59' 
+00°29'48' 
+00°29'42' 
+00°25'57' 
+00° 19^40' 
+00°46'23' 
-00°46'l6' 
+00°3l'08' 
+00° 31^54' 
+00°40'21 
+00°20'30' 
+00° 43' 
+00°42'l9' 
+00° 34^13' 
+00°27'00' 
+00°22'22' 
+00° 34^09' 
+00° 34^09' 
-00°47'20' 
-00°54'l6' 
+00°20'06' 
+00° 19^34' 
+00° 19^35' 
+00° 19^39' 
-00°5l'42' 
+ 00°38''22' 





> 23 


06 




> 22.35 


21.43 ± 0.12 


20.87 ± 0.16 


20.32 ± 0.11 


19 


39 


± 





01 


IS 


46 ± 


18 


46 


18.13 ± 0.01 


17.98 ± 0.01 


17.91 ± 0.01 


18 


OS ± 0.02 


IS 


22 


± 





01 




27 lb 






19 97 i 03 


19 


79 ± 0.03 




18 


SO ± 0.03 








































20 


05 


± 





02 


20 


62 ± 


20 


62 


20.11 ± 0.04 


20.02 ± 0.04 


19.56 ± 0.05 


19 


16 ± 0.04 


IS 


80 


± 





01 




> 22 


37 




21.98 ± 0.20 


22.17 ± 0.26 


21.09 ± 0.22 


20.61 ± 0.14 


19 


76 


± 





03 


19 


52 ± 


19 


52 


18.90 ± 0.02 


18 


62 ± 0.01 


18.18 ± 0.02 


17.83 ± 0.01 


18 


12 


± 





01 




> 21 


79 




> 21.40 




> 20. 


94 


> -99.00 




> 19. 


20 


21 


25 


± 





06 


22 


07 ± 


22 


07 


21.55 ± 0.10 


20.99 ± 0.08 


20.49 ± 0.11 


20.08 ± 0.09 


19 


27 


± 





01 




> 21 


91 




> 21.40 




> 20. 


94 


> 19.61 




> 19. 


20 


20 


85 


± 





04 




> 22 


54 




> 21.73 




> 21. 


05 


> 19.74 




> 19. 


20 


20 


06 


± 





02 






























20 


41 


± 





06 






























20 


56 


± 





07 




> 23 


06 




> 22.35 


19 


70 ± 0.03 


19.38 ± 0.04 


19.09 ± 0.04 














> 22 


43 




> 22.02 


21 


94 ± 0.21 


> 20.99 


20.84 ± 0.17 


19 


95 


± 





02 


19 


06 ± 


19 


06 


18.85 ± 0.02 


18 


62 ± 0.01 


18.39 ± 0.02 


18 


11 ± 0.02 












21 


11 ± 


21 


11 


20.80 ± 0.05 


20.49 ± 0.05 


20.15 ± 0.08 


19 


54 ± 0.06 


18 


87 


± 





01 




> 21 


85 




> 21.40 




> 20. 


94 


> 19.61 




> 19. 


20 


21 


21 


± 





05 


18 


76 ± 


18 


76 


18.39 ± 0.01 


17 


98 


± 


0.01 


17.57 ± 0.01 


17 


22 


± 


0.01 


17 


62 


± 





01 


19 


63 ± 


19 


63 


19.24 ± 0.02 


19 


05 


± 


0.02 


18.64 ± 0.02 


18 


32 


± 


0.02 


18 


52 


± 





01 


21 


17 ± 


21 


17 


20.64 ± 0.06 


20 


19 


± 


0.04 


19.64 ± 0.06 


19 


11 


± 


0.04 


18 


65 


± 





01 




> 22 


37 




> 22.02 


21 


90 


± 


0.20 


> 20.99 


20 


49 


± 


0.12 


19 


29 


± 





02 


19 


81 ± 


19 


81 


19.59 ± 0.03 


19 


28 


± 


0.02 


18.82 ± 0.03 


18 


46 


± 


0.02 


18 


20 


± 





01 


18 


82 ± 


18 


82 


18.58 ± 0.01 


18 


24 


± 


0.01 


17.92 ± 0.01 


17 


67 


± 


0.01 












21 


90 ± 


21 


90 


21.38 ± 0.12 


21 


04 


± 


0.09 


20.94 ± 0.19 


20 


43 


± 


0.11 


19 


22 


± 





01 






























20 


56 


± 





04 


18 


83 ± 


18 


83 


18.59 ± 0.01 


18 


20 


± 


0.01 


17.80 ± 0.01 


17 


53 


± 


0.01 


17 


76 


± 





01 


18 


93 ± 


18 


93 


18.84 ± 0.01 


18 


48 


± 


0.01 


18.20 ± 0.02 


17 


94 


± 


0.01 


18 


22 


± 





01 


19 


83 ± 


19 


83 


19.69 ± 0.03 


19 


33 


± 


0.02 


18.95 ± 0.03 


18 


66 


± 


0.02 












19 


94 ± 


19 


94 


19.69 ± 0.03 


19 


61 


± 


0.03 


19.30 ± 0.04 


18 


98 


± 


0.03 












21 


55 ± 


21 


55 


21.01 ± 0.08 


20 


54 


± 


0.06 


19.82 ± 0.07 


19 


17 


± 


0.04 


18 


69 


± 





01 


20 


47 ± 


20 


47 


20.17 ± 0.04 


19 


86 


± 


0.03 


19.27 ± 0.04 


18 


95 


± 


0.03 


19 


09 


± 





01 


21 


29 ± 


21 


29 


21.16 ± 0.07 


20 


61 


± 


0.06 


20.22 ± 0.09 


19 


73 


± 


0.07 












18 


27 ± 


18 


27 


18.09 ± 0.01 


17 


80 


± 


0.01 


17.52 ± 0.01 


17 


48 


± 


0.01 


17 


43 


± 





01 




> 23 


06 




> 22.35 


20 


60 


± 


0.06 


20.12 ± 0.08 


19 


73 


± 


0.07 














> 23 


06 




> 22.35 


21 


64 


± 


0.16 


> 21.13 


20 


70 


± 


0.15 


19 


80 


± 





03 




> 22 


43 




> 22.02 


22 


11 


± 


0.25 


> 20.99 


20 


91 


± 


0.18 


19 


94 


± 





02 


21 


66 ± 


21 


66 


21.42 ± 0.12 


21 


41 


± 


0.13 


> 20.99 


20 


65 


± 


0.14 


19 


73 


± 





01 






























19 


48 


± 





03 






























19 


96 


± 





02 


20 


82 ± 


20 


82 


20.58 ± 0.04 


20 


50 


± 


0.05 


20.25 ± 0.09 


19 


99 


± 


0.08 


19 


91 


± 





02 


22 


65 ± 


22 


65 


22.23 ± 0.26 


21 


50 


± 


0.16 


20.84 ± 0.18 


20 


28 


± 


0.10 












21 


98 ± 


21 


98 


21.85 ± 0.18 


21 


56 


± 


0.17 


20.99 ± 0.20 


20 


74 


± 


0.16 


20 


36 


± 





04 




> 21 


86 




> 21.40 




> 20. 


94 


> 19.61 




> 19. 


20 


21 


32 


± 





10 






























21 


36 


± 





08 


21 


32 ± 


21 


32 


21.16 ± 0.09 


20 


56 ± 0.06 


20.36 ± 0.11 


19.84 ± 0.07 


19 


27 


± 





01 






























20 


27 


± 





03 


20 


22 ± 


20 


22 


19.94 ± 0.04 


19 


55 ± 0.03 


19.08 ± 0.04 


18 


67 


± 0.02 


18 


73 


± 





01 






























20 


21 


± 





03 



19.04 ± 0.01 
18.28 ± 0.01 



18. SI ± 0.01 

19.65 ± 0.01 

20.83 ± 0.04 
19.21 ± 0.01 
20.39 ± 0.03 

20.21 ± 0.03 
20.08 ± 0.02 

18.66 ± 0.01 
19.63 ± 0.01 

18.22 ± 0.01 
19.20 ± 0.01 
19.44 ± 0.03 

17.68 ± 0.01 

18.94 ± 0.01 



18.20 ± 0.01 
19.15 ± 0.01 
20.33 ± 0.03 
17.92 ± 0.01 



18.99 
19.24 
19.05 
19.04 
19.56 
17.76 
18.09 
19.48 
19.79 
19.93 
19.47 
19.73 
19.87 
19.11 
20.83 



± 0.01 
± 0.01 
± 0.01 
± 0.01 
± 0.04 
± 0.01 
± 0.01 
± 0.01 
± 0.02 
± 0.02 
± 0.04 
± 0.02 
± 0.02 
± 0.01 
± 0.07 



20.72 ± 0.10 
19.63 ± 0.04 
20.12 ± 0.05 
18.95 ± 0.01 
19.95 ± 0.03 



47 


80 


± 


6 


50 


41 


33 


± 


7 


27 


14.29 ± 9.27 





992 





993 








3 


35 


42 


07 


± 


6 


55 


24 


56 


± 


7 


29 


10.37 ± 9.30 





927 





949 





38 


± 0.07 


5 


75 


46 


36 


± 


6 


45 


27 


64 


± 


7 


15 


30.30 ± 9.12 





998 









39 


± 0.07 


1 


44 


46 


26 


± 


6 


46 


34 


17 


± 


7 


25 


20.42 ± 9.32 





989 












1 


92 


46 


74 


± 


6 


56 


30 


97 


± 


7 


16 


6.70 ± S.99 





354 





399 








4 


45 


46 


74 


± 


6 


66 


30 


97 


± 


7 


16 


6.70 ± S. 99 





643 





598 








1 


50 


45 


96 


± 


6 


66 


23 


21 


± 


7 


15 


2.63 ± 9.01 





995 









29 


± 0.08 


4 


67 


45 


76 


± 


6 


61 


38 


65 


± 


7 


14 


32.42 ± 8.92 





896 





977 








1 


85 


45 


12 


± 


6 


36 


37 


32 


± 


7 


15 


15.07 ± S.93 





953 





966 











51 


46 


08 


± 


6 


52 


28 


66 


± 


7 


30 


29.26 ± 9.02 





985 





990 








1 


14 


44 


96 


± 


6 


53 


38 


06 


± 


7 


18 


25.49 ± 8.90 





836 












2 


61 


45 


21 


± 


6 


42 


39 


54 


± 


7 


23 


18.88 ± 8.97 





613 





478 








2 


57 


45 


21 


± 


6 


42 


39 


54 


± 


7 


23 


18.88 ± 8.97 





378 





515 








2 


40 


39 


85 


± 


6 


38 


21 


33 


± 


7 


15 


11.52 ± 9.10 









998 





16 


± 0.17 


2 


24 


44 


81 


± 


6 


48 


44 


94 


± 


7 


15 


34.21 ± 9.12 





882 





881 








6 


31 


44 


02 


± 


6 


53 


21 


20 


± 


7 


40 


9.62 ± 8.99 









928 





25 


± 0.03 


3 


57 


42 


74 


± 


6 


48 


17 


48 


± 


7 


14 


14.55 ± 8.99 





922 





919 





68 


± 0.12 


3 


57 


43 


97 


± 


6 


43 


48 


32 


± 


7 


15 


25.88 ± 8.83 









891 








3 


12 


42 


94 


± 


6 


51 


61 


39 


± 


7 


18 


46.14 ± 8.76 





918 















57 


44 


18 


± 


6 


41 


19 


63 


± 


7 


23 


-0.42 ± 9.33 





959 





948 





25 


± 0.03 


5 


73 


44 


71 


± 


6 


44 


18 


68 


± 


7 


27 


9.49 ± 9.12 





984 









35 


± 0.03 


5 


16 


40 


89 


± 


6 


65 


30 


78 


± 


7 


31 


14.69 ± 8.98 





998 





998 





63 


± 0.10 


2 


13 


42 


99 


± 


6 


59 


43 


26 


± 


7 


16 


15.88 ± 8.91 





987 












4 


12 


38 


04 


± 


6 


61 


28 


16 


± 


7 


14 


5.02 ± S.95 





961 









51 


± 0.12 


2 


68 


42 


06 


± 


6 


40 


11 


64 


± 


7 


16 


9.00 ± 9.00 









975 





11 


± 0.02 


6 


28 


39 


87 


± 


6 


41 


45 


74 


± 


7 


28 


31.52 ± 8.86 





996 





996 








2 


10 


42 


01 


± 


6 


55 


32 


10 


± 


7 


21 


20.21 ± 8.92 





932 





985 








2 


56 


41 


50 


± 


6 


41 


19 


09 


± 


7 


30 


-12.39 ± 9.13 





997 





996 





12 


± 0.02 


2 


33 


41 


25 


± 


6 


53 


22 


19 


± 


7 


17 


1.67 ± 9.10 





967 









21 


± 0.03 


2 


15 


41 


37 


± 


6 


53 


20 


20 


± 


7 


25 


9.76 ± 8.79 









976 





29 


± 0.07 


5 


75 


41 


18 




6 


55 


34 


62 


± 


7 


14 


20.07 ± 9.18 









963 





25 


± 0.06 


5 


53 


41 


09 


± 


6 


50 


20 


53 


± 


7 


20 


-1.89 ± 8.79 





983 





975 





69 


± 0.08 


2 


59 


39 


63 


± 


6 


46 


15 


19 


± 


7 


24 


-3.99 ± 8.95 





882 





944 





63 


± 0.08 


3 


62 


39 


95 


± 


6 


38 


28 


10 


± 


7 


IS 


10.01 ± 9.26 









984 





46 


± 0.18 


3 


02 


3S 


19 


± 


6 


48 


12 


71 


± 


7 


23 


5.37 ± 9.19 





999 





999 





09 


± 0.00 


3 


02 


40 


61 


± 


6 


49 


36 


92 


± 


7 


22 


23.08 ± 8.92 









963 








4 


03 


37 


23 


± 


6 


49 


44 


93 


± 


7 


14 


32.16 ± 8.77 





974 





986 








2 


22 


39 


10 


± 


6 


41 


27 


54 


± 


7 


22 


27.09 ± 9.20 





987 





987 








3 


67 


38 


39 


± 


6 


40 


30 


19 


± 


7 


IS 


20.33 ± 8.92 





890 





889 








6 


24 


37 


71 


± 


6 


61 


26 


04 


± 


7 


17 


24.95 ± 9.08 





983 





986 








2 


79 


38 


79 


± 


6 


55 


31 


27 


± 


7 


24 


14.04 ± 8.90 





545 





545 











64 


38 


79 


± 


6 


66 


31 


27 


± 


7 


24 


14.04 ± 8.90 





452 





452 





46 


± 0.15 


1 


60 


37 


63 


± 


6 


62 


32 


91 


± 


7 


19 


29.28 ± 9.11 









996 








2 


37 


39 


29 


± 


6 


62 


20 


79 


± 


7 


25 


4.80 ± 9.07 





764 





590 








6 


51 


37 


82 


± 


6 


46 


33 


80 


± 


7 


11 


19.53 ± 9.10 





892 












1 


98 


34 


85 


± 


6 


33 


45 


90 


± 


7 


21 


32.27 ± 8.92 





075 





326 








2 


35 


34 


85 


± 


6 


33 


45 


90 


± 


7 


21 


32.27 ± 8.92 





357 





219 





33 


± 0.15 




15 


34 


85 


± 


6 


33 


45 


90 


± 


7 


21 


32.27 ± 8.92 





557 





445 








3 


19 


37 


35 


± 


6 


50 


12 


10 


± 


7 


18 


8.14 ± 9.29 





996 





997 





41 


± 0.04 


3 


12 


38 


05 


± 


6 


62 


23 


97 


± 


7 


16 


3.68 ± 9.22 





950 





926 








1 


25 



Note. - * ID with color-magnitude metliod. " Coordinates are the position of IRAC counteiparts. Photometry from the VIKING 
(Sutherland et al. in prep). The reported SPIRE 350 and 500 fluxes are measured at the postions of 250 ^m sources; no SNR cut is applied. 
4.5 fim from the LR method (section[3Jj. Redshifts are from the SDSS catalogue for sources detected in that survey. ^ The separation between 



survey (Z, Y, J, H, Kb) is extracted in 2" diameter apertures 
R is reliability of the counterparts calculated separately at 3.6 and 
the centroid of the SPIRE emission and the IRAC counterpart. 



TABLE 2 - continued. 



lAUID 


H-ATLAS 
ID 




(J2000) 




(J2000) 






(mag) 






y 

(mag) 






(mag) 




(mag 






(mag) 




(mag) 






4.5 ^im 
(mag) 






•S250'^ 
(mJy) 






•5350'' 
(mJy) 




S5OO 
(mJy 


c 




(3.6,^m) 


(4.5 /jm) 










Scpiiralion ^ 
(arcsec) 


J090855. 2+003051 


4739 




^08^55^ 


08 


+ 00° 30 


50' 


' 7 


20 


16 ± 


20 


16 


19 


88 


± 





03 


19 


69 ± 


0.03 


19 


29 


± 





04 


18.96 ± 0.03 


18 


96 


± 





01 


19 


19 


-i- 





01 


37 


35 


+ 


6 


56 


17 


17 


+ 


7 


17 


6.70 ± 


9. 


03 


Q 


997 


Q 


995 





27 







06 


2 


60 


J090842.2+001815 


4779 


09^ 


^08^42- 


51 


+ 00° 18 


18' 


'2 


19 


72 + 


19 


72 


19 


66 


± 





03 


19 


36 ± 


0.02 


19 


01 


± 





03 


18.69 ± 0.02 


18 


84 


± 





01 












37 


31 


+ 


6 


54 


20 


36 


+ 


7 


18 


15.89 ± 


9 


04 


Q 


953 









36 


± 





07 


5 


10 


J09083 1.3+002251 


4849 


09^ 


^08^31^ 


31 


+ 00° 22 




'2 


21 






08 


20 


85 


± 




07 


20 






20 


04 


± 





08 


19.51 ± 0.05 


19 


04 


± 





01 


19 


39 


-[- 


^ 


03 


36 




-)- 




53 


30 




-[- 




21 


10.68 ± 


S 


78 








983 





51 


± 











J090833.9+004616 


4850 


09^ 


^08^33^ 


93 


+ 00° 46 


16' 


'8 


21 


31 ± 


21 


31 


21 


01 


± 





08 


20 


93 ± 


0.08 


20 


16 


± 





09 


19.66 ± 0.06 












18 


60 


+ 


Q 


01 


36 


38 


+ 


6 


55 


33 


06 


+ 


7 


31 


4.89 ± 


9. 


28 






Q 


999 















90 


J090847. 1+003247 


4930 


09" 


^08^47' 


07 


+ 00° 32 




's 


17 


50 ± 


17 


50 


17 


34 


± 





01 


17 


03 ± 


0.00 


16 


75 


± 





01 


16.54 ± 0.00 


16 


53 


± 





01 


16 


33 


+ 


Q 


01 


37 


94 


± 




53 


22 


44 


+ 


7 


24 


15.05 ± 


8 


92 


Q 


998 









16 


± 





02 


1 


36 


J09093 1.4+003936 


5050 


09" 


"^09^31^ 


53 


+ 00° 39 


38' 


'1 








































> 23. 


00 


20 


07 


+ 


Q 


05 


37 


04 






55 


41 


24 


+ 


7 


23 


28.78 ± 


8 


89 






Q 


989 












1 


61 


J09074 1.5+002 103 


5082 


09" 


^07^41^ 


34 


+ 00° 2 1 




'8 


21 


07 ± 


21 


07 


20 


64 


± 





06 


20 


39 ± 


05 


19 


96 







08 


19.64 ± 0.06 


19 




-j- 
















36 




-j- 




48 


23 




-[- 




30 


21.40 ± 


9 


24 




973 









68 


± 





15 






J090757.3+003006 


5157* 


09" 


"^07^57^ 


73 


+ 00° 30 


05' 


'0 


21 


44 ± 


21 


44 


21 


01 


± 





08 


20 


94 ± 


09 


20 


41 


± 





12 


20.15 ± 0.09 


19 










19 


63 


-[- 


^ 


02 


35 


18 


-j- 




51 


40 


62 


-[- 


ij! 


21 


30.04 ± 


9 


12 




462 




277 





52 







12 


^ 


43 


J090846.2+004110 


5258 


09" 


^08^46' 


38 


+ 00° 41 


ll' 


'0 


20 


03 ± 


20 


03 


19 


67 


± 





03 


19 


28 ± 


0.02 


18 


81 







03 


18.40 ± 0.02 


18 


57 


-)- 




01 


IS 


51 


+ 


Q 


01 


35 


60 


-]- 




57 


17 


51 


+ 


7 


16 


-4.48 ± 


9 


. 17 


Q 


998 


Q 


998 





19 


± 





03 


1 


79 


J09093 8.9+003240 


5267 


09" 


^09^39^ 


00 


+ 00° 32 




'5 


20 






07 


19 


82 


± 





02 


19 


58 ± 




19 


15 


± 





03 


18.89 ± 0.03 


18 


61 


-j- 




01 


18 




-[- 




01 


35 




-)- 






18 




-[- 




24 


16.48 ± 


9 


06 













59 


± 











J090840.6+001725 


5327 


09" 


^08^40^ 


92 


+00° 17 


24' 


'8 


19 


27 ± 


19 


27 


19 


19 


± 





02 


18 




01 


18 


54 


± 





02 


18.29 ± 0.02 


18 


58 


-)- 




01 












34 


91 


-)- 




27 


19 


10 


-[- 




14 


2.43 ± 


9. 


07 




935 









20 


± 





03 




56 


J090759.7+002033 


5366 


09" 


■^07^59' 


67 


+ 00° 20 


34' 


'9 


21 


73 ± 


21 


73 


21 


39 


± 





12 


21 


12 ± 


10 


20 


86 


± 





18 


20.32 ± 0.10 


19 


49 






01 












36 


61 


-j- 




37 


19 


92 






26 


2.40 ± 


9. 


03 




996 


















49 


J090932.9+002242 


5386 


09" 


■^09^33^ 


01 


+ 00° 22 


45' 


'6 






































21 


38 


i 




09 


20 


54 


-[- 




05 


34 


12 


-|- 




43 


39 


10 


-[- 




26 


27.32 ± 


8 


96 




824 




958 














23 


J0903 14.0+004235 


5422 


09" 


^03^ 13^ 


98 


+ 00° 42 


3l' 


'2 


19 


61 ± 


19 


61 


19 


39 


± 





02 


19 


-^0 ± 


02 


19 


00 


± 





03 


18.78 ± 0.03 


18 


97 


-]- 




01 


18 


94 


-[- 




01 


33 


19 


-j- 




55 


32 


95 


-[- 




19 


4.46 ± 


8. 


98 




981 




995 





24 


± 





04 




85 


J090914.8+002041 


5450 


09" 


"^09^14^ 


82 


+ 00° 20 


39' 


'5 


21 


62 ± 


21 


62 


21 


43 


± 





09 


21 


07 ± 


0.09 


20 


62 


± 





13 


20.47 ± 0.13 


20 


46 


-)- 




03 


20 


85 


+ 


Q 


05 


35 


94 


+ 


6 


42 


24 


72 


+ 


7 


23 


16.42 ± 


9 


06 





987 





922 





39 


± 





19 


2 


44 


J090929.6+003313 


5521 


09" 


^09^29^ 


72 


+ 00° 33 




'2 












> 21. 


73 












> 19. 


74 


> 19.13 


19 


64 


-j- 




01 


19 




-[- 




01 


35 




-)- 




31 


24 




-[- 




15 


8.91 ± 


8. 


99 
























J090837.1+005002 


5530 


09" 


^08^37^ 


04 


+ 00° 50 


02' 


'2 


20 


91 ± 


20 


91 


20 


67 







06 


20 


10 ± 


0.04 


19 


72 







06 


19.21 ± 0.04 












18 


17 


+ 


Q 


01 


36 


26 


± 


g 


56 


16 


39 


+ 


7 


27 


— 0.12 ± 


9 


.02 









999 





60 


± 





06 





93 


J090825.8+004217 


5538 


09" 


"^08^25^ 


83 


+ 00° 42 




'7 


22 






80 


22 


16 


± 





23 


21 


76 ± 




21 


21 


± 





24 


20.68 ±0.14 


20 




^ 






20 




-[- 




03 


35 








41 


17 




-[- 




32 


— 5.24 ± 


9 


.06 




736 




















J09082 1.9+002445 


5564 


09" 


^08^22- 


14 


+ 00° 24 


47' 


'3 




> 22 


37 




22 


06 


± 





21 


21 


45 i 


0.13 


20 


77 


± 





16 


20.13 ± 0.09 


19 


30 






01 


19 


18 


± 





01 


35 


10 


± 




62 


19 


57 


+ 


7 


16 


7.32 ± 


8. 


97 





992 





993 












3 


29 


J090910.9+003517 


5621 


09" 


"^09^ 11^ 


35 


+ 00° 35 




'4 


20 


81 ± 




81 


20 


71 







05 


20 


34 ± 




20 


13 


± 





08 


19.66 ± 0.06 


18 


97 






01 


19 




-J- 




01 


33 




± 




56 


25 




-[- 




15 


26.89 ± 


8 


96 













71 


± 





11 






J090848. 1+002626 


5691 


09" 


^08^48^ 


13 


+ 00° 26 


22' 


' 1 


20 


10 ± 


20 


10 


19 


89 


± 





03 


19 


57 ± 


0.03 


19 


19 


± 





04 


18.86 ± 0.03 


19 


04 


-[- 




01 


19 


29 


-|- 





01 


35 


24 


± 




38 


18 


50 


+ 


„ 


16 


16.50 ± 


8 


90 





885 





890 





39 


± 





05 


4 


15 


J090812.1+002430 


5735 


09^ 


^08^ 12^ 


26 


+ 00° 24 




'7 


19 


21 ± 




21 


18 


78 


± 





01 


18 






18 


45 


± 





02 


18. 16 ± 0.02 


18 


31 


+ 





01 


18 




-[- 




01 


34 




-)- 




51 


15 




-[- 




25 


3. 13 ± 


8. 


99 













16 













J090852.9+004106 


5774 


09" 


^08^52' 


90 




07' 


'1 






85 






> 21. 


40 




>'20 


9^4 




> 19. 


61 




> 19.20 


19 


82 


± 





02 


19 


58 


-[- 




01 


33 


83 


-)- 




45 


38 


19 


-[- 




17 


22.59 ± 


9 


26 




996 




996 














79 


J090912.3+002129 


6012 


09^ 


^09^ 12^ 


65 


+ 00° 2 1 


28' 


'0 


19 


82^± 


19 


82 


19 


59 







02 


19 


27 ± 


02 


18 


90 


± 





03 


18.63 ± 0.03 


18 


38 


± 





01 


18 


81 


± 




01 


35 


13 


-j- 




65 


24 


35 


-[- 




22 


6.36 ± 


9. 


07 




997 




995 





56 


± 





02 




93 


J090908.3+002545 


6100* 


09^ 


^09^08^ 


64 


+ 00° 25 


48' 


'4 


22 


45 i 


22 


45 


21 


84 


± 





13 


21 


55 zt 


0.14 


21 


16 


± 





21 


20.48 ± 0. 13 


19 


52 


± 





01 


19 


35 




Q 


01 


33 


08 


± 




38 


44 


67 


+ 


7 


17 


42.48 zt 


8 


86 





244 





540 












4 


78 


J090902.6+004737 


6146 


09^ 


^09^02^ 


42 


+ 00° 47 




'5 


21 






63 


21 


56 


± 





13 


20 






20 


48 


± 





12 


19.89 ± 0.07 












19 




-|- 




01 


32 




-j- 




37 


20 




-[- 




18 


2.26 ± 


8. 


98 













66 







22 






J090928. 7+002630 


6189* 


09^ 


^09^29^ 


03 


+ 00° 26 


31' 


'0 




> 22 


54 






> 21. 


73 






> 21 . 


05 




> 19. 


74 


> 19.13 


20 


90 


± 





05 


20 


48 


+ 


Q 


03 


33 


47 


+ 


6 


52 


35 


08 


+ 


7 


28 


35.56 ± 


8 


97 


Q 


302 


Q 


430 












4 


06 


J090754.6+003345 


6220* 


09" 


^07^54' 


82 


+ 00° 33 


42{ 


'2 


19 


71 ± 


19 


71 


19 


41 


± 





02 


19 


14 ± 


0.02 


18 


79 


± 





03 


18.48 ± 0.02 


IS 


68 


± 





01 


18 


77 







01 


34 


05 


± 


6 


51 


12 


37 




7 


23 


8.23 ± 


9. 


14 





518 





570 





19 


± 





04 


4 


19 


J090946.4+003847 


6224 


09" 


^09^46^ 


39 


+00°38 


47' 


'1 


20 


48 ± 


20 


48 


20 


19 


± 





03 


19 


83 ± 


0.03 


19 


25 


± 





04 


18.81 ± 0.03 


18 


65 


± 





02 


18 


76 


± 





01 


34 


77 


± 


6 


59 


14 


14 


± 


7 


24 


4.09 ± 


9. 


06 





998 





999 





46 


± 





12 





97 


J090832.0+002749 


6324 


09" 


'^08^31- 


94 


+00°27 


52^ 


'4 


22 


04 ± 


22 


04 


21 


85 







17 


21 


59 ± 


0.15 


20 


77 


± 





16 


20.22 ± 0.10 


19 


18 


± 





01 


19 


21 


± 





01 


33 


98 


± 


6 


35 


23 


50 


± 


7 


19 


10.73 ± 


9 


12 





989 





972 












2 


90 


J090956.0+003739 


6348 


09" 


^09^56- 


03 


+00°37 


37^ 


'2 


19 


46 ± 


19 


46 


19 


18 







02 


18 


77 ± 


0.01 


18 


31 


± 





02 


17.86 ± 0.01 


18 


22 


± 





01 


18 


IS 


± 





01 


33 


90 




6 


45 


12 


66 


± 


7 


24 


-3.81 ± 


9 


.38 





999 





999 





28 


± 





06 


2 


12 


J090840. 1+001928 


6480 


09" 


^08^40^ 


02 


+00° 19 


30' 


'1 


19 


88 ± 


19 


88 


19 


80 







03 


19 


50 ± 


0.02 


19 


37 


± 





05 


19.03 ± 0.03 


19 


25 


± 





01 












33 


32 


± 


6 


55 


25 


13 


± 


7 


23 


27.36 ± 


9 


19 





995 









31 


± 





09 


2 


65 


J090847.1+001512 


6792 


09" 


^08^47^ 


28 


+ 00° 15 


14' 


'1 


18 


37 ± 


18 


37 


18 


16 


± 





01 


17 


79 ± 


0.01 


17 


45 


± 





01 


17.14 ± 0.01 


17 


55 


± 





01 












33 


60 


± 


6 


38 


6. 


58 




7. 


14 


-1.24 ± 


8 


.87 





999 









26 


± 





05 


2 


71 


J0903 24.0+003954 


6893 


09" 


^03^23^ 


95 


+00°39 


52{ 


'1 


21 


76 ± 


21 


76 


21 


35 


± 





11 


20 


98 ± 


0.09 


20 


48 


± 





13 


19.86 ± 0.07 


19 


17 


± 





01 


19 


32 


± 





01 


32 


93 


± 


6 


49 


16 


27 


± 


7 


32 


-0.91 ± 


9 


.29 





946 





955 





49 


± 





19 


2 


54 


J090808.2+002115 


6962 


09" 


"^08^08' 


14 


+00°21 


16' 


'5 


22 


38 ± 


22 


38 


21 


56 


± 





13 


21 


07 ± 


0.10 


21 


30 


± 





27 


20.14 ± 0.09 


19 


12 


± 





01 












32 


57 


± 


6 


48 


17 


35 


± 


7 


23 


-4.52 ± 


9 


.05 





997 
















1 


28 


J09 1308.7-005605 


6991 


09" 


^081^08- 


80 


-00°56 


08^ 


'0 




> 22 


58 






> 21. 


33 






> 20. 


91 




> 19. 


72 




> 19.23 


20 


36 







03 


20 


03 







02 


30 


39 


± 


6 


44 


38 


53 




7 


14 


27.37 zt 


8 


94 





987 





986 












2 


47 


J090915.4+003450 


7107 


09" 


^09^15^ 


06 


+ 00° 34 


48' 


'8 


20 


68 ± 


20 


68 


20 


31 


± 





04 


19 


99 ± 


0.03 


19 


41 


± 





04 


19.00 ± 0.03 


18 


62 







01 


18 


81 







01 


32 


15 




6 


41 


8. 


28 




7. 


23 


1.37 ± 


9. 


01 





972 





976 





74 


± 





11 


5 


75 


J090855. 7+0021 21 


7216 


09" 


^08^55^ 


64 


+00°21 


23^ 


'1 


21 


87 ± 


21 


87 


21 


44 


± 





12 


21 


04 ± 


0.09 


20 


82 


± 





17 


20.38 ± 0.11 


19 


34 


± 





01 


18 


88 


± 





01 


32 


34 


± 


6 


39 


37 


62 


± 


7 


16 


33.56 zt 


8 


89 





996 





998 












2 


02 


J090836.5+002513 


7826 


09" 


"^08^36^ 


54 


+00°25 


09' 


'5 


22 


11 ± 


22 


11 


21 


54 


± 





13 


21 


13 ± 


0.10 


21 


16 


± 





23 


20.33 ± 0.11 


19 


41 


± 





01 


19 


28 


± 





01 


30 


94 


± 


6 


33 


35 


84 


± 


7 


13 


27.94 ± 


9 


06 





991 





992 












3 


62 


J090842.8+003623 


8022 


09" 


^08^43^ 


03 


+00°36 


2l{ 


'7 




> 21 


85 






> 21. 


40 




> 20. 


94 




> 19. 


61 




> 19.20 


20 


96 


± 





06 


20 


45 


± 





03 


29 


66 


± 


6 


41 


56 


85 


± 


7 


21 


50.43 ± 


9 


18 





970 





981 












3 


07 


J091309.5-005417 


9262 


09" 


"^08^09^ 


85 


-00°54 


19' 


'3 






































20 


42 


± 





02 


20 


04 


± 





02 


28 


12 


± 


6 


57 


39 


70 


± 


7 


23 


19.64 ± 


9 


22 





958 





954 












4 


46 


J090830.0+002439 


9290 


09" 


^08^29^ 


79 


+00° 24 


38^ 


' 1 




> 21 


85 






> 21. 


40 




> 20. 


94 




> 19. 


61 




> 19.20 


19 


93 


± 





02 


19 


72 


± 





01 


29 


27 


± 


6 


47 


38 


29 


± 


7 


26 


19.79 ± 


8 


98 





754 





770 












3 


45 


J090909.4+003440 


10110 


09" 


"^09^09^ 


80 


+ 00° 34 


38' 


'5 




> 22 


54 






> 21. 


73 






> 21. 


05 




> 19. 


74 


> 19.13 


21 


29 







05 


20 


62 


± 





04 


27 


62 


± 


6 


52 


37 


27 


± 


7 


23 


27.12 ± 


9 


02 





558 





861 












5 


06 


J090820.3+002808 


11918 


09" 


^08^19^ 


84 


+00°28 


08' 


'6 


22 


34 ± 


22 


34 


21 


49 


± 





13 


21 


16 ± 


0.10 


20 


69 


± 





15 


20.53 ± 0.13 


19 


95 


± 





03 


19 


99 







03 


26 


08 




6 


45 


37 


76 


± 


7 


28 


16.30 zt 


9 


11 





764 





787 












6 


92 


J090834. 7+003635 


15947 


09" 


^08^34^ 


77 


+00°36 


39' 


' \ 






































20 


55 


± 





03 


20 


27 


± 





02 


22 


52 


± 


6 


50 


37 


86 


± 


7 


14 


32.28 ± 


9 


03 





731 





852 












3 


43 



Note. - * ID with color-magnitude method. " Coordinates are the position of IRAC counterparts. * Photometry from the VIKING survey {Z,Y, J, H, Ks) is extracted in 2" diameter apertures 
(Sutherland et al. in prep). The reported SPIRE 350 and 500 fim fluxes are measured at the posticus of 250 fim sources; no SNR cut is apphed. R is reliability of the counterparts calculated separately at 3.6 and 
4.5 fim from the LR method (section im . Redshifts are from the SDSS catalogue for sources detected in that survey. ^ The separation between the centroid of the SPIRE emission and the IRAC counterpart. 



